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pattern  shape  of  electromagnetic  radiation  from  flush  mounted  antennas 
used,  in  high  speed  aircraft,  missiles  or  satellites. 
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I.  INTRODUCTION 


This  report  deals  with  a geometrical  theory  of  diffraction  (GTD)  [1] 
analysis  of  the  high  frequency  (h.f.)  radiation  from  slots  in  perfectly- 
conducting,  planar  and  circular  cylindrical  structures  which  are  partly 
covered  with  an  impedance  surface  patch.  The  present  study  is  applicable 
to  the  problem  of  controlling  the  pattern  shape  of  the  electromagnetic 
radiation  from  a flush  mounted,  airborne  antenna  for  satellite  coitmuni- 
cation  purposes.  In  this  application,  a slot  antenna  is  flush  mounted  in 
the  aircraft  fuselage,  and  the  radiation  pattern  in  the  roll  plane 
containing  the  slot  is  of  interest.  Consequently,  the  problem  is 
essentially  two-dimensional  (2-D)  in  nature,  and  the  aircraft  fuselage  may 
be  modeled  approximately  by  a circular  cylinder  in  the  roll  plane.  The 
effects  of  the  aircraft  wings  are  ignored  in  this  initial  study;  however, 
their  effects  can  be  incorporated  in  the  future  work  in  a straightforward 
manner  as  will  be  indicated  later  on.  In  particular,  it  is  of  prime 
interest  in  the  satellite  communication  type  application  for  the  airborne 
antenna  to  have  a gain  approximately  that  of  the  pattern  maximum  near  the 
horizon  or  the  shadow  boundary  associated  with  the  antenna.  Such  a desired 
high  gain  may  be  achieved  by  an  impedance  loading  of  the  perfectly- 
conducting  surface  around  the  slot. 

The  antenna  configurations  of  interest  which  are  analyzed  in  this 
report  are  illustrated  in  Figures  1 and  2.  The  slot  is  present  in  the 
perfectly-conducting  surface  and  the  perfectly-conducting  surface  is 
coated  with  a sufficiently  long  impedance  surface  patch  which  also  covers 
the  slot  as  shown  in  Figures  1 and  2.  With  a proper  choice  of  impedance, 
the  slot  may  be  allowed  to  radiate  as  an  end-fire  antenna  near  the  horizon; 
this  effect  may  then  be  employed  to  increase  the  gain  near  the  horizon. 

By  a proper  choice  of  the  impedance,  Zj,  it  is  implied  that  Zs  is  chosen 
such  that  the  slot  can  excite  a bound  surface  wave  mode  on  the  planar 
impedance  surface  of  Figure  1;  whereas,  it  can  excite  an  Elliot  type  mode 
[2]  on  the  curved  impedance  surface  of  Figure  2.  The  Elliot  mode  propa- 
gates along  the  curved  impedance  surface  with  negligible  leakage,  and  it 
corresponds  in  the  limiting  case  of  the  infinite  radius  of  curvature  to 
the  bound  surface  wave  mode  on  the  planar  impedance  surface.  These  modes 
on  the  impedance  surface  diffract  from  the  ends  of  the  impedance  patch 
(as  a result  of  the  discontinuity  in  surface  impedance  there),  thereby 
producing  an  end-fire  effect.  It  is  noted  that  the  surface  wave  or 
Elliot  type  modes  discussed  above  are  excited  by  an  axial  (or  z-directed) 
slot  in  the  configurations  of  Figures  1 and  2,^if  the  surface  impedance 
is  inductive.  Whereas,  a circumferential  (or  i-directed)  slot  will  excite 
these  modes  only  if  the  surface  impedance  is  capacitive. 

In  the  present  analysis,  the  impedance  surface  patch  is  used  to 
approximately  simulate  the  effects  of  a thin,  uniform  dielectric  or 
ferrite  cover,  or  a properly  designed  corrugated  surface  of  finite 
extent  on  the  perfectly-conducting  planar  or  cylindrical  antenna 
structure.  The  slot  in  the  perfect! y conducting  structure  is  covered 
by  the  dielectric  or  ferrite  material,  and  it  is  assumed  that  the 
electric  field  in  the  slot  aperture  is  known,  so  that  one  also  knows 
the  equivalent  magnetic  current  in  the  aperture.  The  radiation 
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Figure  1.  A thin,  2-0  slot  in  an  infinitely  long,  planar, 
perfectly-conducting  structure  with  a surface  ' 
impedance  patch  of  length  L.  The  value  of  the 
surface  impedance  is  Zg. 
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Picture  2.  A thin  2-'!  l .L  in  a perfectly-conducti ng  circular 
cylin  ^.,  rtith  an  impedance  surface  patch  of  arc 
length  L.  The  value  of  the  surface  impedance  is  Zg. 
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from  a slot  of  finite  width  may  be  obtained  by  quantizing  the 
source  distribution  in  the  aperture,  and  by  superposing  the  fields 
radiated  by  each  of  the  quantized  sources.  An  equivalent  quantized 
source  in  this  case  is  a z-directed  magnetic  line  source  for  the  axial 
slot;  whereas,  it  is  a T-directed  magnetic  line  dipole  for  the  circum- 
ferential slot.  These  equivalent  quantized  magnetic  currents  radiate  in 
the  presence  of  the  perfectly-conducting  surface  (which  covers  the 
aperture  as  well),  and  the  dielectric  or  ferrite  material  which  covers 
the  perfectly-conducting  surface.  The  strengths  of  the  equivalent 
sources  are  weighted  according  to  the  source  distribution  in  the  aperture. 
In  the  present  report,  only  the  radiation  from  a magnetic  line  source  and 
a magnetic  line  dipole  will  be  analyzed  since  the  radiation  from  a slot 
of  finite  width  may  be  readily  obtained  from  this  analysis  by  super- 
position. In  the  impedance  surface  approximation  for  the  truncated 
dielectric  or  ferrite  cover  on  the  perfectly-conducting  surface,  the 
equivalent  problem  to  be  analyzed  consists  of  an  equivalent  magnetic  line 
or  line  dipole  source  on  the  impedance  surface  patch  which  covers  the 
perfectly-conducting  planar  or  cylindrical  surface. 

While  the  cylindrical  antenna  configuration  of  Figure  2 is  important 
in  this  study  because  it  models  an  aircraft  fuselage  in  the  roll  plane, 
the  planar  configuration  of  Figure  1 is  also  analyzed  in  this  report 
because  it  is  a more  basic  geometrical  configuration  involving  a 
perfectly-conducting  structure  with  a surface  impedance  patch,  and  because 
a comparison  of  the  patterns  for  the  configurations  in  Figures  1 and  2 
allows  one  to  ascertain  the  effects  of  surface  curvature  especially  in 
the  region  of  space  in  which  the  slot  is  directly  visible  (i.e.,  the  lit 
region) . 

A GTD  analysis  has  been  performed  previously  with  success  to 
analyze  the  radiation  from  slots  in  truncated,  planar  dielectric  covered 
surfaces  [3],  hence  it  was  decided  to  extend  the  GTD  analysis  of  [3]  to 
treat  the  problems  in  Figures  1 and  2.  In  particular,  the  GTD  solutions 
to  the  problems  in  Figures  1 and  2 are  "built  up"  from  the  asymptotic  h.f. 
solutions  to  several,  appropriate  canonical  problems  as  in  [3],  however 
the  type  of  canonical  problems  analyzed  and  employed  in  this  report  are 
quite  different  from  the  ones  employed  in  [3].  The  significant  advantages 
of  the  GTD  method,  which  is  an  asymptotic  h.f.  method,  are  that  once  the 
pertinent  GTD  launching,  propagation  and  diffraction  coefficients  are 
known  from  the  solutions  to  canonical  problems  (or  by  some  other  means), 
then  the  GTD  may  be  used  to  calculate  the  fields  radiated,  scattered,  or 
diffracted  from  complex  structures  in  a simple  fashion  by  appropriately 
expressing  the  fields  in  terms  of  the  above  mentioned  coefficients 
together  with  the  spatial  divergence  factors  for  the  rays  associated  with 
these  fields.  Such  a ray  description  results  from  the  local  nature  of 
propagation,  scattering,  and  diffraction  which  is  exhibited  by  asymptotic 
h.f.  field  approximations.  The  spatial  divergence  factors  for  the  rays 
indicate  the  presence  and  the  locations  of  the  caustics  of  the  ray 
system.  General Iv  on„  oi  the  caustics  of  the  diffracted  ray  occurs  at 
the  point  of  diffraction  on  the  surface.  Away  from  the  point  of  dif- 
fraction, the  diffracted  ray  tube  spreads  according  to  the  laws  of 
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ordinary  geometrical  optics.  The  launching  or  diffraction  coefficients 
provide  the  starting  amplitudes  of  the  rays  directly  emanated  from  the 
source,  or  the  rays  diffracted  from  certain  points  of  diffraction  on  the 
surface,  respectively.  Therefore,  in  addition  to  the  usual  rays  of 
geometrical  optics  corresponding  to  the  incident  and  reflected  rays  which 
propagate  from  the  source  to  the  observation  point  along  direct  and 
reflected  ray  paths,  respectively,  the  diffracted  rays  propagate  from  the 
source  to  the  observation  point  via  points  of  diffraction  on  the  surface 
in  accordance  with  the  generalized  Fermat’s  principle  proposed  in  Keller's 
GTD  [1].  The  diffracted  rays  are  produced  by  the  presence  of  discontin- 
uities or  variations  in  the  geometrical  and  electrical  properties  of  the 
surface  (e.g.,  a discontinuity  in  surface  impedance).  For  a comprehensive 
account  of  GTD  which  also  includes  some  new  and  useful  results  for 
engineering  applications,  one  is  referred  to  Kouyoumjian  [4].  The  ray 
description  for  the  fields  not  only  provides  a simple  solution,  but  one 
which  is  also  highly  efficient  for  numerical  computations.  Furthermore, 
the  ray  contributions  provide  a direct  physical  insight  into  the  various 
scattering  diffraction  and  radiation  mechanisms  which  are  present, 
thereby  also  providing  useful  design  information  for  controlling  these 
effects  advantageously. 

One  notes  however  that  the  ray  optical  field  description  breaks 
down  at  and  in  the  neighborhood  of  caustics  and  shadow  boundaries, 
hence  the  ray  field  must  either  be  modified  or  supplemented  by  a 
separate  solution  so  that  it  remains  valid  within  these  regions  of  space 
which  are  commonly  referred  to  as  transition  regions.  The  field  analyses 
in  the  transition  region  is  complicated  because  the  fields  must  change 
rapidly  but  smoothly  across  these  regions.  Modifications  or  the  intro- 
duction of  supplemental  solutions  to  the  ray  field  within  the  transition 
regions  are  referred  to  as  uniform  representations  if  they  not  only 
remain  valid  within  these  transition  regions,  but  if  they  also  reduce  to 
the  proper  ray  fields  outside  the  transition  regions.  There  are  no  ray 
caustics  in  the  far  zone  of  the  sources  in  the  antenna  configurations  of 
Figures  1 and  2,  respectively;  however,  shadow  boundaries  are  present  for 
the  circular  cylinder  antenna  configuration  of  Figure  2 due  to  the 
finite  radius  of  curvature  of  the  cylinder  in  this  case.  Such  shadow 
boundaries  are  absent  in  the  planar  configuration  of  Figure  1.  Hence,  a 
uniform  asymptotic  approximation  is  developed  in  the  present  analysis 
for  the  cylindrical  antenna  configuration  of  Figure  2 to  provide  field 
expressions  which  are  valid  within  the  transition  regions,  and  which  reduce 
uniformly  to  the  ray  fields  outside  these  regions. 

The  analytical  details  are  presented  in  Section  II.  Some  preliminary 
numerical  results  for  the  radiation  patterns  are  illustrated  in  Section 
III  for  both,  axial  and  circumferential  slots  in  the  geometrical  configur- 
ations of  Figures  1 and  2,  respectively.  The  solutions  to  the  various 
canonical  problems  from  which  the  final  GTD  solution  is  constructed  are 
indicated  in  the  AppendiC".  Some  useful  modifications,  refinements, 
and  generalization:  of  this  work  which  are  worth  pursuing  are  also 
discussed  in  Section  III. 
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II.  ANALYTICAL  FORMULATION 


The  ray  optical  analysis  based  on  the  GTD  for  the  antenna  configur- 
ations of  Figures  1 and  Z will  be  indicated  in  this  section.  Specifically, 
the  manner  in  which  the  analysis  of  the  composite  problems  in  Figures  1 
and  2 will  be  built  up  or  constructed  from  the  asymptotic  solutions  to 
somewhat  simpler  canonical  problems  will  be  indicated.  Basically,  the 
total  radiated  field  is  composed  of  the  fields  associated  with  rays  which 
propagate  not  only  directly  from  the  source  to  the  observation  point,  but 
also  with  r-ays  which  propagate  via  points  on  the  surface  at  which 
variations  or  discontinuities  in  the  geometrical  and  electrical  properties 
can  occur.  Such  points  constitute  points  of  diffraction  on  the  surface. 

The  planar  configuration  of  Figure  1 will  be  analyzed  first  in 
Section  A,  and  the  analysis  for  the  cylindrical  configuration  of  Figure 
2 will  follow  in  Section  B.  Both,  the  axial  and  circumferential  slot 
cases  will  be  considered  in  Sections  A and  B.  It  is  noted  that  due  to 
the  2-0  nature  of  the  problem,  an  axial  slot  in  the  configurations  of 
Figures  1 and  2 radiates  only  a z component  of  the  magnetic  field  which 
will  be  denoted  by  whereas,  a circumferential  slot  in  these  configur- 
ations radiates  only  a z component  of  the  electric  field  which  will  be 
denoted  by  Ez.  The  solution  for  the  axial  slot  case  may  be  referred  to 
as  the  TEz  solution,  and  likewise  the  solution  for  the  ci ri_umferential 
slot  case  may  be  referred  to  as  the  TMz  solution. 

As  mentioned  previously  in  Section  I,  only  the  radiation  from  a 
magnetic  line  source  (for  the  TE^  case)  and  the  magnetic  line  dipole 
source  (for  the  TMz  case)  on  an  impedance  surface  patch  which  covers  a 
perfectly-conducting  planar  or  circular  cylindrical  surface  will  be 
analyzed  in  this  report;  these  magnetic  line  currents  constitute  the 
quantized,  equivalent  sources  in  the  axial  and  circumferential  slot 
apertures.  The  radiation  from  a 2-D  axial  or  circumferential  slot  of 
finite  width  may  be  readily  obtained  from  this  analysis  via  superposition. 
In  the  impedance  surface  approximation  for  the  truncated  dielectric  or 
ferrite  cover  on  the  perfectly-conducting  surface,  the  equivalent 
problems  to  be  analyzed  therefore  consist  of  an  equivalent  magnetic  line 
or  line  dipole  source  on  an  impedance  surface  patch  which  covers  the 
perfectly-conducting  planar  and  cylindrical  surfaces,  respectively. 

An  eJ“^  time  convention  is  assumed  and  suppressed  in  the  following 
analysis  unless  specified  otherwise  (as  in  Appendix  II). 

A.  Analysis  of  the  Planar  Antenna  Configuration 

According  to  the  GTD,  one  may  simply  describe  the  field  radiated 
by  the  planar  antenna  configuration  of  Figure  1 in  terms  of  a super- 
position of  the  fields  associated  with  the  rays  which  are  directly 
radiated  by  the  source  at  0,  and  the  rays  which  are  diffracted  from  the 
ends  of  the  surface  impedance  patch  at  Qi  and  Q2,  respectively.  The 
pertinent  rays  for  this  problem  are  illustrated  in  Figure  3 in  which 
u*"  represents  the  field  of  the  ray  dirpctlv  radiated  from  a line  source 
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AT  0 
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Figure  3.  Rays  associated  with  the  radiation  from  a line  source 
on  an  infinite,  planar  perfectly-conducting  surface 
which  is  covered  with  an  impedance  surface  patch. 


at  0,  whereas  u-|  and  U2  represent  the  fields  diffracted  from  the  ends 
Ql  and  Q2  of  the  surface  impedance  patch,  respectively.  The  surface 
impedance,  Zs  is  chosen  to  support  a bound  surface  wave  type  mode  which 
is  xcited  by  the  line  source  at  0.  Thus  the  surface  impedance,  Zs  is 
inductive  for  the  magnetic  line  source  case  (corresponding  to  the  axial 
slot,  or  TE2  case),  whereas  Zs  is  capacitive  for  the  magnetic  line 
dipole  source  (corresponding  to  the  circumferential  slot,. or  TM^  case). 

The  field  of  this  bound  surface  wave  mode  is  denoted  by  ui ; this  field 
propagates  from  the  source  at  0 to  the  ends  of  the  surface  impedance 
patch  at  Ql  and  Q2,  respectively  as  shown  in  Figure  3,  Since  there  is  a 
discontinuity  in  surface  impedance  at  Qi  and  Q2,  the  field  yi  at  Qi ^2  is 
partly  reflected,  and  partly  diffracted  from  Qi,2-  Thus,  u^  2 is  produced 

via  the  diffraction  of  the  surface  wave  field  ui(Qi  o)  when  it  impinges 
upon  the  surface  impedance  discontinuity  at  Qi  2-  course,  multiple 
reflections  of  the  surface  wave  field  can  occur  between  the  ends  Qi  and  Q2, 
and  the  inclusion  of  these  multiple  refl ection-diffraction  effects  is 
discussed  later  in  Section  III.  It  is  assumed  in  the  present  analysis 
that  the  surface  impedance  patch  is  sufficiently  large  in  extent  so  that 
only  the  bound  surface  wave  mode  field  u'  (which  is  a non-radiating  field) 
constitutes  the  dominant  contributor  to  the  field  on  the  impedance 
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surface  far  from  the  source,  since  the  field  u*"  (i.e.,  the  radiated  or 
space  wave  field)  becomes  vanishingly  small  along  the  surface  impedance 
boundary  when  it  is  evaluated  far  from  the  source.  This  assumption  may  be 
verified  by  solving  the  canonical  problem  of  a line  source  on  an  infinite, 
planar  impedance  surface  as  shown  in  Figure  4.  The  canonical  problem 
of  Figure  4 is  briefly  analyzed  in  Appendix  I,  and  the  solution  to. this 
problem  provides  the  necessary  expressions  for  the  fields  u*"  and  u^  of 
Figure  3. 


P 


Figure  4,  Geometry  of  the  canonical  problem  of  a line  source  on 
an  infinite,  planar  impedance  surface. 


In  the  radiation  pattern  analysis,  the  field  point  P is  assumed 
to  be  in  the  far  zone  of  the  surface  impedance  patch.  The  field  u*"(P), 
which  is  the  field  u*"  evaluated  at  P is  a cylindrical  wave  type  field 
which  emanates  from  the  line  source  at  0 when  P is  sufficiently  far  from 
0;  this  field  may  be  represented  according  to  geometrical  optics  as 


u’'^(P)  -V  u’'(P') 


\|  p'"+[PP'] 


,-jk[PP'] 


(1) 


where  P'  is  some  reference  point  on  the  ray  path  from  0 to  P,  and  k is 
the  free  space  wave  number.  The  caustic  distance  of  this  ray  is 
simply  the  distance  OP'.  In  order  to  relate  the  field  u>"(P)  to  the 
source,  one  must  take  the  limit  of  the  RHS  of  Equation  (1)  as  p>"-K)  (or 
P'  approaches  the  point  0).  In  the  ray  approximation,  u'"(P')  is  singular 


quantity  PP'  in  Equation  (1)  denotes  the  distance  from  P' 


to  P. 
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at  the  caustic,  but  the  limit  of  u'^(P')  as  p'^-^0  exists,  and  it  is 
proportional  to  the  strength  of  the  magnetic  line  source  or  the  magnetic 
line  dipole  source  at  0 and  the  radiation  coefficient  R which  is  defined 
below.  Let  R = uM  6(x)6(y)  denote  this  magnetic  line  source,  or  the 
magnetic  line  dipole  source  so  that 


lim  ^ u^(P')  a M R . (2) 


One  notes  that  u=z  for  the  axial  slot  or  the  magnetic  line  source  case, 
whereas  u=x  for  the  circumferential  slot  or  the  magnetic  line  dipole  case. 
The  constant  of  proportional i ty  in  the  above  limit  is  defined  as  Cq,  where 
Co  is  some  complex  constant.  The  radiation  coefficient  R indicates  the 
manner  in  which  the  source  M at  0 distributes  the  radiation  field  in  space 
thus,  R is  identical  to  the  pattern  factor  of  the  source.  In  general,  R 
is  a function  of  the  azimuthal  angle  4>.  One  notes  that  PP'->p  as  p’"-K),  so 


u'"(P)  ^ M R(^) 


(3) 


The  constant  Cq  and  the  radiation  coefficient  R((p)  are  determined  from 
the  solution  to  the  canonical  problem  of  Figure  4,  as  given  in 
Appendix  I. 


for  the  magnetic  line  dipole,  or  TM^  case 
for  the  magnetic  line  source,  or  TE^  case.  (4) 


The  Yq  ’n  Equation  (4)  is  the  free  space  admittance.  R(<t>)  will  be  dif- 
ferent "'or  the  magnetic  line  source,  and  the  magnetic  line  dipole;  the 
explicit  form  of  this  coefficient  will  be  indicated  later. 


The  fields  u-|(P)  and  U2(P)  at  P which  are  diffracted  from  the  ends 
of  the  surface  impedance  patch  at  Q]  and  Q^,  respectively  may  be 
assumed  to  be  produced  by  localized  equivalent  magnetic  line  sources  or 
line  dipoles  at  these  ends.  Thus,  u<J  o(P)  ^Iso  be  expressed  in 
terms  of  a cylindrical  ray  divergence ’factor  via  ray  optics  as 


(5) 


I I 

where  P-i  is  some  reference  point  between  Q]  and  P;  likewise  P2  is  some 
reference  point  between  Q2  and  P.  The  caustic  distances  p]  2 these 
diffracted  ray  fields  are  simply  given  by  p^j*  2 = Q]  2 2-'  before. 
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is  finite  and  proportional  to  the  strength  of  the  equivalent  line  source 
at  Qi  2 The  strength  of  the  equivalent  line  source  would  of  course  be 
proportional  to  the  strength  of  the  surface  wave  field  u’’  which  strikes 
the  ends  Qi  ^2  to  produce  the  diffracted  field  u^  tience 


1 im 


ri.2 


u^Qi^2^ 


D (41 1 2 ) » 


where  D(4-|  p)  is  the  diffraction  coefficient  at  Q-|  2 which  indicates  the 
manner  in  Anich  the  diffracted  field  is  distributed  in  space.  The 
functional  form  of  D is  the  same  at  Q]  and  Q2,  respectively,  and  the 
diffraction  coefficient  at  Q]  depends  on  the  angle  4i;  whereas,  the  dif- 
fraction coefficient  at  Q2  depends  on  the  angle  42-  The  angles  4]  and  42 
are  indicated  in  Figure  3,  Thus,  Equation  (5)  becomes 


-jks 


1.2 


1,2 


= 71-4  . 


(6) 


d d 

It  is  noted  that  P]P"^s^  as  p°-*0,  and  P2^'*S2  Equation  (b).  In 

the  far  zone,  one  may  approximately  replace 


-jks 


1.2 


by 


’1.2 


~7=~  ■ ® 

> P 


using  the  origin  as  the  phase  reference.  Thus, 


jk^os4i^2  e ^ ^ 


(7) 


► p 


The  incident  surface  wave  field  u^(0]  launched  by  the  source  M at  0 
may  be  represented  as 


u’(Qi,2)  = M e 


-i4 


(B) 
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sw 

where  L and  g are  the  surface  wave  launching  and  propagation  coef- 
ficients, respectively;  both  of  these  quantities  are  found  from  the 
solution  to  the  canonical  problem  of  Figure  4 as  given  in  Appendix  I. 
The  remaining  quantity,  D(i|ii  o)  Equation  (7)  must  be  found  from  the 
solution  to  the  canonical  problem  of  Figure  5 as  given  in  Appendix  II. 


Z,#0 


z,=o 


u*  = INCIDENT  SURFACE  WAVE  MODE  FIELD 
u'’>  REFLECTED  SURFACE  WAVE  MODE  FIELD 
DIFFRACTED  FIELD 


Figure  5.  Geometry  of  the  canonical  problem  of  surface  wave 
diffraction  by  a di sconti nui ty  in  surface  imped- 
ance associated  with  a two  part,  planar  impedance 
surface  in  which  the  surface  impedance  for  x>0 
and  y=0  is  zero  corresponding  to  a perfectly- 
conducting  boundary. 


As  pointed  out  by  Keller  and  Karal  (Journal  of  Appl . Phys.,  Vol.  31, 

No.  6,  June  1960,  pp.  1039-1046),  and  more  recently  by  Felsen  and 
Choudhary  (IEEE  Trans.  Antennas  and  Propagation,  Vol.  AP-21 , November 
1973,  pp.  827-842),  the  geometrical  representation  of  evanescent  or 
surface  waves  is  by  complex  rays;  however,  since  the  source  and  field 
points  lie  on  the  lossless  impedance  surface  in  the  calculation  of 
‘^MQ1,2)  this  analysis,  and  since  the  field  of  the  surface  waves  is 
not  included  at  P,  the  representation  in  Equation  (8)  is  adequate. 

One  may  now  specialize  the  above  set  of  results  to  the  axial  slot 
(or  TEz)  and  the  circumferential  slot  (or  Tf^)  cases,  respectively.  As 
mentioned  previously,  the  magnetic  current  M corresponding  to  a 
quantized  equivalent  source  in  the  slot  is  given  by 


z M 6(x)6(y)  for  the  magnetic  line  source  corresponding 

to  the  axial  slot  or  TMz  case. 

X M 6(x)6(y)  for  the  magnetic  line  dipole  corresponding 

to  the  ci rcumferential  slot  or  TE^  case. 
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in  which  M is  the  amplitude  of  the  source  which  is  assumed  known  and 
'5(x)6(y)  is  the  two-dimensional  Dirac  delta  function.  The  corresponding 
fields  for  the  two  slot  types  are 


and  0 < < TT. 


00) 


The  source  R in  Equation  (9)  for  the  TE^  and  TM^  cases  is  explicitly 
defined.  The  quantity  R((|>)  in  u*^  for  the  TE^  and  TM2  cases  is  given 
separately  in  Equations  (A-26)  and  (A-45),  respectively  of  Appendix  I. 

The  diffraction  coefficient  D((ji)  appearing  in  uy  2 above  is  presented 
separately  for  the  TEz  and  TMz  cases  in  Equation^  (A-74)  and  (A-97), 
respectively  of  Appendix  II.  One  recalls  that  u^  2 ^^so  requires  an 
explicit  knowledge  of  u^Ql  2)  which  in  turn  requires  a knowledge  of 
Lsw  and  3 as  mentioned  in  Equation  (8).  These  quantities  are  explicitly 
given  in  Appendix  I by  Equations  (A-27)  and  (A-19)  for  the  TEz  case,  and 
by  Equations  (A-46)  and  (A-40)  for  the  TMz  case. 

B.  Analysis  of  the  Circular  Cylinder  Antenna  Configuration 

There  are  sow  essential  differences  between  the  radiation  mechanisms 
present  in  the  planar  antenna  configuration  of  Figure  1 which  was 
analyzed  in  part  A,  and  the  cylindrical  antenna  configuration  of  Figure  2 
which  is  analyzed  below.  These  differences  will  become  evident  in  the 
following  discussion.  As  mentioned  earlier,  the  magnetic  current  R 
corresponding  to  a quantized  equivalent  source  in  the  slot  aperture  is 
given  by 


M = G M ; 

P 


The  boundary  SBQ'  in  Figure  6 constitutes  the  shadow  boundary  for  the 
source  R at  Q'.  The  field  radiated  by  the  source  M at  Q'  may  be 
obtained  via  GTD  in  terms  of  a s^erposition  of  the  fields  of  the  rays 
directly  radiated  by  the  source  M,  and  the  fields  of  the  rays  which  are 
diffracted  from  the  ends  of  the  surface  impedance  patch  at  Q]  and  Q2, 
respectively  as  shown  in  F;^qure  6^.  The  field  directly  radiated  by  M 
is  denoted  by  whereas  u^  and  u^  represent  the  fields  diffracted  from 
Q-|  and  Q2,  along  direct^ray  paths  to  the  observation  point  P,  as  in 
Figure  6a.  The  fields  u<j*  and  u^  may  be  assumed  to  be  produced  by 
equivalent  magnetic  line  sources  or  line  dipoles  at  Q-)  and  Q2,  respectively. 
Hence,  the  boundaries  SBQ]  and  SBQ2  illustrated  in  Figure  6 correspond 
to  the  shadow  boundaries  of  equivalent  sources  of  the  diffracted  fields 
at  Q]  and  Q2,  respectively. 


z,  for  the  magnetic  line  source 


u = 


or  the  TE^  case. 


(11) 


<t>,  for  the  magnetic  line  dipole 


or  the  TM  case. 
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Figure  6.  Rays  and  shadow  boundaries  associated  with  the 
cylindrical  antenna  configuration  of  Figure  2. 
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As  in  the  planar  case,  the  type  of  impedance  surface  chosen  is  such 
that  it  supports  an  Elliot  type  surface  wave  mode  which  propagates  along 
the  impedance  surface  from  the  slot  to  the  ends  Qi  and  Q2  of  the  patch 
with  negligible  leakage  (provided  th^. cylinder  is  sufficiently  large). 

The  Elliot  mode  field  is  denoted  by  u) , and  it's  propagation  path  is 
illustrated  in  Figure  6.  The  field  u^  impinges  upon  the  termination  of 
the  sur^^ce  impedance  patch  at  Qi  and  Q2  to  produce  the  diffracted 
fields  DfV  and  lW,  respectively  which  traverse  the  direct  ray  paths 
QlP  ^nd  Q2P  to  the  observation  point  P.  However,  in  addition  to  M 
and  u^,  there  exist  an  additional  class  of  diffracted  fields  resulting 
from  the  diffraction  of  creeping  waves  that  are  excited  at  Qi  and  Q2; 
these  creeping  waves  propagate  over  the  perfectly-conducting  portion 
of  the  cylinder.  Sufficiently  far  from  Q]  and  Q2  the  creeping  waves 
are  expressed  in  terms  of  a rapidly  converging  set  of  surface  ray  modes 
whose  amplitudes  decay  exponentially  due  to  a continuous  shedding  of  rays 
from  the  surface  ray  mode  fields  (via  surface  diffraction)  along  the 
forward  tangents  to  their  surface  ray  paths  on  the  perfectly-conducting 
part  of  the  cylinder.  These  surface  ray  modes  were  introduced  by 
Keller  [1]  to  analyze  the  diffraction  by  smooth,  convex  cylinders  via 
GTD.  The  surface  ray  modes  are  also  commonly  referred  to  as  Wa;^son 
modes  [5].  The  fields  of  the  surface  ray  modes  are  denoted  by  u<^w 

whereas,  the  fields  of  the  rays  diffracted  from  these  surface  ray 
mod^s  are  denoted  by  u^“  and  2f|° , respectively.  The  rays  corresponding 
to  u^^  and  u^'"'  travel  in  opposite  directions  along  the  perfectly- 
conducting  portion  of  ;^e  cylinder  as  shown  in  Figure  6.  The  diffracted 
rays  corresponding  to  u^°  and  u|“  are  also  illustrated  in  Figure  6.  In 
Mrticular,  one  notes  that  the  rays  corresponding  to  the  diffracted  field 
uy^2  eJ^ist  on  the  lit  side  of  the  shadjjw  boundary  SBQi^2’  whereas,  the  rays 
corresponding  to  the  diffracted  field  u?“p  exist  on  the  shadow  side  of 

SBQ1.2. 

The  ray  descriptions  U]  2 ^nd  u]^2  on  either  side  of  the  shadow 
boundary  SBQi  ^2  and  in  the  neighborhood  of  SBQi  2-  Therefore,  a 

uniform  approximation  for  the  field  which  remains  valid ’within  this 
;^ansition  region  adjacent  to  SBQi^2i  and  which  reduces  to  u'f  2 and 
Ui°2  outside  the  transition  region’in  the  lit  and  shadow  side^  of 
SB0i,2»  respectively  is  presented  in  terms  of  the  Fock  integrals  [5]. 

One  notes  that  in  addition  to  the  Elliot  mode  field  excited  by  TT 
on  the  curved  impedance  patch,  there  also  exist  a set  of  Watson  type 
modes  which  leak  off  energy  via  diffraction_from  the  curved  impedance 
boundary  as  they  propagate  from  the  source  M at  Q'  to  the  ends  of  the 
impedance  boundary  at  Q]  and  Q2.  The  existence  of  these  Watson  type 
modes  is  not  indicated  in  Figure  6 because  their  diffraction  effects 
will  be  assumed  to  be  negligible  in  comparison  with  those  due  to  the 
Elliot  type  mode.  The  Elliot  mode  propagates  with  very  little  leakage, 
hence  it  may  be  assumed  that  this  mode  is  the  dominant  contributor  to 
the  field  on  the  curved  impedance  surface  which  is  incident  at  Qi,2* 

While  the  assumption  that  the  leakage  from  the  Watson  type  modes  Seing 
small  may  be  true  only  for  a certain  range  of  parameters  (such  as  those 
involving  tne  value  of  the  surface  impedance,  Zs ; the  length  of  the 
curved  impedance  surface  section;  and  the  curvature  of  the  cylinder). 
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the  present  analysis  is  based  on  this  assumption  being  true*.  There  is 


also  a transition  region  adjacent  to  the  shadow  boundary  SBQ'  associated 
with  the  source  M at  Q*.  One  notes  that  the  ray  field  u’^  far  from  Q' 
vanishes  along  SBQ';  this  may  be  verified  via  an  asymptotic  solution  to 
the  canonical  problem  of  Figure  7 as  given  in  Appendix  III.  Thus,  since 


X 


Figure  7.  Geometry  of  the  canonical  problem  of  a line  source 
(or  thin  slot)  on  a circular  cylinder  which  is 
completely  covered  with  an  impedance  surface. 


the  leakage  or  surface  diffraction  effects  of  the  Watson  type  modes  is  I 

neglected,  and  u>^  is  continuous  everywhere  (even  near  the  boundary  SBQ'  I 

where  it  vanishes)  on  the  lit  side  of  SBQ'  (u*^  0 on  the  dark  or  shadow  J 

side  of  SBQ'),  no  uniform  approximation  for  the  fields  is  presented 
within  the  transition  region  adjacent  to  SBQ'.  However,  one  would  have 
to  modify  u*^  if  the  leakage  of  the  Watson  modes  from  the  impedance  surface 
is  included.  This  modification  would  provide  a uniform  approximation  with- 
in the  transition  region  adjacent  to  SBQ';  such  a uniform  approximation 
involves  a modified  Fock  integral.  One  notes  that  the  inclusion  of  the 
Watson  modes  also  requires  that  their  diffraction  from  the  ends  of  the 
impedance  patch  at  Q]  and  Q2  be  included.  Furthermore,  a uniform  ap- 
proximation to  the  diffraction  of  these  Watson  modes  from  Qi^2  '’’ust  then 
be  included  within  the  transition  regions  adjacent  to  SBQi  2^  otherwise, 

a discontinuity  in  the  field  pattern  would  result  across  s5Qi^2-  I 

inclusion  of  the  effects  of  these  Watson  type  modes  on  the  curved  impedance  I 

section  is  proposed  as  part  of  the  future  work  dealing  with  refinements  I 

and  extensions  to  the  present  work. 

By  dividing  the  space  surrounding  the  cylinder  in^o  six  regions  as 
indicated  in  Figure  6,  one  notes  that  the  total  field  u radiated  by  the 
slot  may  be  given  within  each  of  these  regions  I,  II,  III,  IV,  V,  and  VI, 
respectively,  in  terms  of  the  ray  fields  u’",  u“,  and  u^“2  follows. 

2 

*0ne  migi;t  expect  this  assumption  to  be  fairly  accurate  for  sufficiently 

large  cylinders  and  for  impedance  patch  lengths  which  are  small  in 

comparison  to  the  size  of  the  cylinder.  , 


14 


a 


(12) 


I 

\ 

t 


As  noted  previously,  u = in  the  TM^  problem;  whereas, u = in  the 
TE2  problem, 

^CW  ^cw 

The  fields  u-|  and  up  impinge  upon  the  discontinuities  in  surface 
impedance  at  Qz  and  Q-j,  after  being  launched  from  Q-)  and  Qp,  respectively, 
thereby  giving  rise  to  further  reflec;^ions  and  diffractions.  However, 
these  diffraction  and  reflections  of  u'f'^  and  from  Qp  and  Q] , 
respectively  are  neglected  in  the  present  analysis  since  the  fields 
uj”p  are  exponentially  decaying  fields,  and  for  sufficiently  large 
excursions  of  the  ray  paths  from  Qi  to  Qp,  and  vice  versa,  over  the 
perfectly  conducting  portion,  the  fields  if^'^p  are  small  enough  to  where 
their  diffraction  effects  at  Q2,l  become  eveA  smaller.  The  cylinder  is 
taken  to  be  sufficiently  large  electrically  in  the  present  analysis, 
so  that  the  diffraction  of  Dfy'^p  from  Qp  1 is  indeed  negligible  in  com- 
parison to  the  other  diffraction  effects  being  considered.  Also,  multiple 
reflections  of  the  Elliot  mode  field  £(''  can  occur  between  Qi  and  Qp  on 
the  impedance  surface.  The  inclusion  of  these  multiple  reflection- 
diffraction  effects  is  discussed  later  in  Section  III. 

The  fields  u*^  and  u^  whose  ray  paths  are  indicated  in  Figure  6 may 
be  found  from  the  solution  to  the  canonical  problem  of  Figure  7 as  given 
in  Appendix  III.  Following  the  representation  for  Of*"  for  the  flat  surface 
case  of  part  A,  one  may  likewise  express  ray-optical ly  as 


u^P)  M RU) 


P in  lit  region  of  Q' . 


(13) 


The  distance  s'  = Q'P,  and  Cq  has  been  defined  earlier  in  Equation  (4). 
The  radiation  coefficient  R(ip)  in  Equation  (13)  is  the  same  as  that  in 
Equation  (3)  for  the  planar  case;  this  result  could  have  been  directly 
predicted  via  the  geometrical  optics  ray  approximation  without  having 
to  solve  for  £(’^  as  in  Appendix  III.  Choosing  the  origin  as  the  phase 
reference  allows  one  to  express  Of'"  as 


15 


(14) 


il’^(P)  'V  M R(<p)  — 


where  the  angle  ip  is  illustrated  in  Figure  6. 


Next,  the  Elliot  mode  field  u^  on  the  curved  impedance  surface 
which  arrives  at  Q]  or  Q2  from  Q'  may  also  be  expressed  in  a fashion 
similar  to  ul  for  the  planar  case  (see  Equation  (8))  as 


^ (^1.2) 


M L®”  e 


t^  ^2  = length  Q'Q^  ^ 


0W 

The  quantities  L and  y are  the  Elliot  mode  launching  and  propagation 
coefficients,  respecti vely , which  are  found  in  Appendix  III. 

The  field  u^  which  is  incident  at  Q]  2 produces  diffracted  rays 
which  emanate  from  Q-i^2  as  mentioned  earlier.  For  a sufficiently  large 
cylinder,  the  nature  of  the  diffraction  of  the  Elliot  mode  from  Qi ^2 
may  be  assumed  to  be  locally  the  same  as  that  in  Figure  5 for  the  planar, 
two-part  impedance  surface,  except  in  the  transition  region  adjacent 
to  the  shadow  boundary  SBQi  2 ^nd  on  the  shadow  side  of  SBQ]^2  below  the 
transition  region  (wherein  {he  surface  diffracted  field  the  only 

ray  field  present).  However,  even  within  the  transition  region  adjacent 
to  SBQ-]  2 and  on  the  shadow  side  of  this  transition  region,  one  may 
indirec{ly  employ  the  surface  wave  diffraction  coefficient  for  the  planar 
case  of  Figure  5 to  obtain  an  equivalent  magnetic  current  source  at  Qi  2 
which  in  turn  produces  u^°2  and  the  transition  region  field.  The 
equivalent  current  is  a magnetic  line  source  for  the  TEz  case;  whereas, 
it  is  a magnetic  line  dipole  for  the  TMz  case.  This  "equivalent"  source 
is  inhomogeneous  since  it  is  a function  of  6-|  or  62,  i.e.,  it  depends  on 
the  aspect,  41,  and  it  is  regarded  as  the  source  which  produces  the 
diffracted  fields.  Let  the  equivalent  magnetic  current  at  Qi  2 be 
designated  as  R®'!,  where  ’ 


pn  <5(p-a)S(4-<('_) 
“1 

PQ  5(p-a)'5(4i+<f  ) 
U2  1^2  (62)  ^ 


; at  Q, 


; at  Qn 


z,  for  TE^  case 
, for  TM^  case 


z,  for  TE  case 

} , (17a. b) 


t2«  for  TM^  case 


.jut'  rjv*" 


in  which  ti  2 '’^  vector  at  Q]^2  shown  in  Figure  6.  The 

strength  z)  equivalent  source  is  evaluated  in 

Appendix  IV’whicfi  deals  with  the  solution  to  the  canonical  problem 
of  Figure  8.  Since  the  details  are  presented  in  Appendix  IV,  only  the 


TO  P ON 

SHADOW  SIDE  OF  SSO, 


Figure  8.  Geometry  of  the  canonical  problem  of  the  diffraction  of 
an  Elliot  mode  by  a discontinuity  in  surface  impedance 
on  a circular  cylinder. 


final  risults  for  the  diffracted  ray  fields  il^  2 and  2fi  o which  are 
produced  by  at  Qi  2 are  indicated  .below.  One  recalls  that  irf  £ exists 
on  the  lit  side, of  SB$i  21  whereas,  Lr^o  exists  on  the  shadow  side’of  SBqi^2- 
The  ray  field  2 spreads  .cylindricalTy  outward  from  Qi  2-  Foi"  the  ^ke 
of  brevity,  the  forms  of  and  will  be  indicated  below,  however  J2 
have  exactly  similar  forms  (and  may  be  obtained  by  replacing  the 
subscript  1 by  2). 


d^(P)  ^ 


P in  the  lit  side  of  SBQ.|  . 

(18) 
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The  distances  si  and  S2  correspond  to  QiP  and  Q2P,  respectively,  as 
in  Figure  6,  ana  the  angles  6]  and  62  are  also  shown  in  Figure  C.  The 
quantity  R®*^(6i)  is  the  radiation  coefficient  of  the  source 


p 


at  Q] , and  it  is  specified  in  Appendix  IV.  However,  since  u-|(P)  is  a 
diffracted  ray  field,  it  is  more  appropriate  to  represent  it  in  terms 
of  a diffraction  coefficient  rather  than  a radiation  coefficient.  In 
fact,  it  is  shown  in  Appendix  IV  that 


C M 


eq  peq. 


p)  = u'(Qi  p)D(6,  ,) 


'o"'l,2  " '"1,2'  ■ “ "^1,2'^^"1,2 

where  D is  same  as  the  diffraction  coefficient  in  Equation  (7).  Thus, 


(19) 


-jks^ 

!rf(P)  'V  lI’(Q  )D(6  ) S ; P in  the  lit  side  of  SBQ, 

(and  0 < 6^  < tt)  (20) 

as  in  Equation  (6)  for  the  planar  case.  In  the  far  zone,  one  may  write 


i/^(P)  'v  &^{Q^)D(6^) 


jkacos(<|i-<(>^) 


(21) 


One  notes  that  u^  0 for  > tt;  this  is  consistent  with  the  assump- 
tion that  th^  curved  surface  diffraction  effects  of  the  Watson  type  modes 
(excited  by  M®9)  on  the  curved  impedance  patch  are  negligible. 

The  surface  diffracted  ray  field  excited  by  M^92(6-|  £=0)  over 
the  perfectly-conducting  surface  is  givert  by  Pathak  and  l^ouyoumjian  [6] 
as 


-[«P+jk]A, 

e ^ 


-jks, 
e 


when  P is  on  the  shadow  side  of  SBQ^.  (22) 

Only  u|^  is  indicated  above,  and  is  similar  to  The  distance 

1-|  2 ■ length  Qi  2^1,2-  distance  S-|  2 " ^1  2^  where  T]  ^ is  the 
^tnt  of  diffraction’or  shedding  of  the  surface  diffracted  ray  field 
u^^2  Ql,2  is  the  launching  point  of  this  surface  diffracted  ray  field 
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i 

i 


I 

t 


as  indicated  in  figure  6.  Lp,  Dp  and  ap  are  the  launching,  diffraction 
and  attenuation  coefficients  of  the  pth  surface  ray  (or  Watson)  mode; 
these  are  discussed  in  [6]  and  are  also  introduced  in  Appendix  IV.  In 
numerical  calculations,  the  inclusion  of  only  the  first  couple  of  modes 
(p=l,2)  is  sufficient  to  obtain  good  accuracy  provided  P is  outside  the 
transition  region  adjacent  to  SBQi- 

One  may  now  specialize  the  above  results  to  the  axial  slot  (or 
TEz)  and  the  ci rcumferenti al  slot  (or  TMz)  cases,  respectively.  The 
magnetic  current  R at  Q'  corresponding  to  the  two  types  of  slots  is 
defined  in  Equation  (11)  in  which  M is  the  amplitude  of  the  source  which 
is  assumed  known.  The  total  field,  iJ  radiated  by  the  source  R is  obtained 
from  Equation  (12)  after  noting  that 


[tm  I 

for  the  ^ ’ cases  ; and  0 ^ _ 2tt  . (23) 


The  total  field  u in  Equation  (12)  is  composed  of  p and 

vXl  "sd 

The  diffracted  fields  ^ u^  ^ Turn  are  expressed  in  terms 
of  lI^(Qi  2)-  Finally,  one  recalls  that  the  expressions  for  2* 

u^'^2»  d^(0^  2)  ^oiTtain  quantities  R(<)>),  L^'^,  y,  M^^2  (°*"  ^'^)t 

D,  Lp,  On  and  ap.  In  particular,  R(<p)  in  u*^  is  given  by  Equations 
(A-125)  and  (A-140)  in  Appendix  III  for  the  TEz  and  TMz  cases,  respec- 
tively. The  and  y are  given  in  Appendix  III  by  Equations  (A-123) 
and  (A-124)  respectively  for  the  TEz  case;  whereas,  they  are  given 
by  Equations  (A-138)  and  (A-139)  respectively  for  the  TMz  case.  The 
strengtn  of  the  equivalent  current  Mj^2(*^l  2^  given  by  Equations 
(A-145)  and  (A-152)  for  the  TEz  ^nd  Tl^z  cases,  respectively  in 
Appendix  IV.  Also  given  in  Appendix  IV  are  the  expressions  for 
Lp,  Dp  and  ctp  for  the  TEz  and  TMz  cases,  respectively;  specifically, 
one  is  referred  to  Equations  (A-143),  (A-144b),  (A-149),  (A-151), 

^-147),  and  ^A-155)  for  these  expressions.  Also,  the  D appearing  in 
iW  2 and  in  (of  Equations  (A-145)  and  (A-152))  is  given  by  Equations 
(A-74)  and  (A-§7)  for  the  TEz  and  TMz  cases,  respectively  in  Appendix  I.'. 


Finally,  one  notes  that  the  ray  fields  u^  2 ^nd  uf^2  on  the  lit 
and  shadow  sides  of  the  shadow  boundary  SBQi^2’are  not  valid  within 
the  transition  region  adjacent  to  SBQ]^2‘  The  angular  extent  of  the 
transition  region  is  0[(ka)"V3]  radians.  Following  a generalization 
of  the  analysis  in  [6]  which  is  based  on  the  results  of  Ivanov[7], 
one  obtains  the  necessary  .modifications  which  constitute  uniform  approx- 
imations for  the  fields  u^  and  uj  in  tne  transition  region  by  firstly 
replacing  the  expression  in  Equation  (22)  for  with 
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uf  (P)  M®'’ 


g(^i) 

g(Ci) 

-jM,  -^'=1 
e ' 1 for 

/Si 

TE 

z 

TM 

'case,  (M)  1 

' 

1 m , 

V 

'-d  .. 


and  by  replacing  the  expression  in  Equations  (18)  or  (20)  for  u-j  with 


u^(P)  . 


-jks. 


for 


TE 


l™z| 


case.  (25) 


Equation  (24)  for  il]  is  to  be  employed  for  P in  the  shadow  zone  of 
the  transition  region;  whereas.  Equation  (25)  for  is  to  be  employed 
for  P in  the  lit  zone  of  the  transition  region.  The  Fock  functions 
g(-) , and  §((•)  are  defined  in  [6];  they  are  also  defined  and  tabulated 
in  [5]. 


mH-i 


1 A. 

C-|  = ; 5 = -m  cos  (^  - '5])  . ('5-[<  ^ in  transition  region) 

(26a, b) 


and 


/ka\ 
m = K-j) 


1/3 


(27) 


Although  the  above  field  expressions  in  Equations  (24)  and  (25)  are 
associated  with  the  transition  region  adjacent  to  SBQ] , similar  express- 
ions apply  to  the  transition  region  adjacent  to  SBQ2.  It  was  mentioned 
previously  that  inhomogeneous,  i.e.,  it  depends  on  the 

aspect  (P  via  the  artgle  6-|^2  (see  Equation  (19)  for  example);  however, 
it  is  important  to  note  tliat  while  wf^2  Equations  (18),  (l9),  and 
(25)  depend  on  6-|  g;  the  Equations  (22)  and  (24)  is  evaluated 

at  4-|^2  0 3S  indicated  in ’detail  in  Appendix  IV. 
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III.  NUMERICAL  RESULTS  AND  DISCUSSION 


Numerical  results  for  the  radiation  patterns  of  the  antenna  con- 
figurations of  Figures  1 and  2 are  presented  in  this  section;  the 
I patterns  have  been  calculated  via  the  GTD  results  which  are  developed 

I for  these  configurations  in  Section  II.  The  radiation  patterns  for  the 

I axial  slot  or  TE^  case  are  presented  for  the  planar  and  cylindrical  con- 

I figurations  in  Figures  9-15;  whereas,  those  for  the  circumferential  slot 

I or  TMz  case  are  presented  in  Figures  16-20.  These  patterns  are  normalized 

I to  0 dB  at  their  peak  value. 


The  results  of  the  GTD  calculations  are  compared  against 
corresponding  results  which  are  obtained  via  an  independent  calculation 
based  on  a combination  of  the  method  of  moments  (MM)  and  the  GTD;  the 
latter  method  abbreviated  as  MM-GTD  as  applied  to  the  antenna  configurations 
of  Figures  1 and  2 is  discussed  in  a separate  report  [8].  The  MM-GTD  is 
also  an  approximate  method;  however,  it  is  expected  to  be  slightly  more 
accurate  than  the  initial  GTD  analysis  presented  in  this  report.  It  is 
noted  that  in  general,  the  results  obtained  by  the  present  GTD  analysis 
agree  well  with  those  obtained  by  the  I'ti-GTD  analysis  of  [8]. 


In  general,  one  notes  that  increasing  the  size  of  the  impedance  patch 
increases  the  number  of  pattern  ripples  in  the  lit  region  (i.e.,  in 
the  region  of  space  where  the  slot  is  directly  visible);  whereas,  increasing 
the  size  of  the  cylinder  tends  to  decrease  the  level  of  the  radiation 
pattern  in  the  deep  shadow  region  of  the  slot.  Over  the  range  of  imped- 
ance values  considered  in  the  present  calculations,  one  notes  that  in- 
creasing the  value  of  the  surface  impedance  for  the  TE^  case  in  general 
increases  the  size  of  the  pattern  ripple  in  the  lit  region;  whereas,  the 
reverse  appears  to  be  true  for  the  TMz  case. 


The  high  lobe  structure  caused  by  the  use  of  the  surface  impedance 
patch  is  a direct  consequence  of  a poor  impedance  match  at  the  edges  of 
the  patch.  This  is  harmful  for  the  TE^  case  but  it  is  a disaster  for 
the  TMz  excitation.  Tapering  for  impedance  matching  purpose  has  been 
used  in  the  other  report  [8]  but  it  would  not  be  of  any  real  import 
as  it  was  used  in  so  far  as  the  TMz  geometry  is  concerned  since  the 
presence  of  the  zero  surface  impedance  would  always  cause  a strong 
reflection.  Consideration  should  be  given  to  terminating  the  structure 
in  an  impedance  surface  whose  value  is  chosen  so  that  no  surface  wave 
exists.  This  should  force  the  energy  off  the  surface  and  in  the  desired 
direction  with  minimum  back  reflection  as  is  the  case  in  corrugated 
horn  structures. 


The  GTD  method  presented  in  this  report  yields  a simple  and 
efficient  procedure  to  analyze  the  radiation  patterns  of  the  antenna 
configurations  in  Figures  1 and  2;  furthermore,  due  to  the  local  nature 
of  the  GTD,  it  provides  a physical  insight  into  the  radiation  mech- 
anisms involved.  The  numerical  results  indicate  that  this  method 
indeed  works  well.  However,  some  refinements  and  modifications  of 
the  present  analysis  are  required  to  make  the  present  GTD  analysis 
useful  and  more  accurate  under  general  conditions.  For  example,  it 
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would  be  worth  investigating  the  following  topics  in  this  connection.  It 
is  of  interest  to  investigate  whether  a more  accurate  representation  for 
the  fields  on  the  impedance  surface  region  would  improve  the  GTD  analysis 
for  Smaller  size  impedance  patches,  and  even  smaller  size  cylinders.  In 
the  planar  case,  one  notes  that  the  possibility  of  multiple  surface  wave 
reflection  - diffraction  exists,  and  this  interaction  could  be  quite 
important  if  the  surface  wave  reflection  coefficient  at  the  ends  of 
impedance  patch  is  large.  When  interactions  are  important,  it  is  im- 
portant to  have  an  accurate  representation  for  the  surface  field  which 
impinges  on  the  ends  of  the  impedance  patch;  otherwise,  the  errors 
could  accumulate  with  the  inclusion  of  each  higher  order  multiple 
interaction.  Furthermore,  for  the  size  of  the  impedance  patches  considered 
in  the  present  calculations  may  not  be  long  enough  to  assume  that  the 
space  wave  component  of  the  field  of  a slot  in  an  impedance  surface  is 
vanishingly  small  at  the  ends  of  the  patch,  hence  some  errors  could 
result  from  not  employing  a more  accurate  representation  for  the  field  on 
the  impedance  surface  in  this  case.  One  notes  that  only  the  surface  wave 
mode  or  the  Elliot  type  mode  field  is  assumed  to  exist  over  the  planar  or 
curved  impedance  sections,  respectively,  in  this  study  however  this  may  be 
true  only  for  large  size  impedance  sections.  Intimately  tied  in  with  this 
improved  surface  field  calculation  is  the  inclusion  of  the  effects  of  the 
Watson  type  modes  on  the  curved  impedance  section  to  obtain  a better  surface 
field  representation  in  this  case;  however,  since  a calculation  of  the 
propagation  constants  of  these  Watson  type  modes  could  present  some 
numerical  problems  for  certain  values  of  the  surface  impedance,  an  alterna- 
tive surface  field  expression  involving  a modified  Fock  integral  for  the 
curved  impedance  surface  should  be  investigated.  Of  course,  one  would 
have  to  look  at  ways  to  approximate  this  Fock  type  integral  for  the  field 
on  the  curved  impedance  surface  for  ease  of  numerical  computations. 
Furthermore,  the  inclusion  of  not  only  the  Elliot  mode,  but  also  the 
Watson  modes  on  the  curved  impedance  patch  leads  to  a more  complicated 
analysis  for  the  diffraction  of  the  Watson  modes  from  the  ends  oi  the 
impedance  patch;  in  particular  the  analysis  of  this  problem  is  more 
complicated  within  the  transition  regions  adjacent  to  shadow  boundaries 
SBQ]  and  SBQ2  (see  Figure  6),  than  that  presented  here  for  only  the 
Elliot  mode  case.  The  diffraction  effects  of  the  Watson  type  modes 
are  assumed  to  be  negligible  in  the  present  analysis;  the  conditions  under 
which  this  assumption  is  true  needs  to  be  more  carefully  investigated. 

It  would  be  interesting  to  study  the  transition  region  problem 
associated  with  the  diffraction  of  the  Watson  type  modes  over  a curved 
impedance  patch,  as  it  is  the  dominant  effect  within  this  transition 
region  for  curved  impedance  sections  which  do  not  support  an  Elliot 
type  mode  (e.g.,  a capacitive  reactance  surface  in  the  axial  slot  or 
TEz  case).  Some  of  these  problems  will  be  investigated  in  the  future 
phases  of  this  study. 

These  proposed  modifications  are  expected  to  further  improve  the 
accuracy  of  the  GTD  calculations  presented  in  this  section,  in  addition 
to  making  the  GTD  method  work  under  somewhat  more  general  conditions 
which  may  be  encountered  in  practical  antenna  problems  of  the  type  shown 
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in  Figures  1 and  2.  While  the  use  of  the  MM-GTD  method  for  analyzing  these 
problems  [8]  is  far  more  efficient  than  the  conventional  MM  techniques, 
it  is  less  efficient  than  the  analysis  based  on  just  the  GTD  alone.  At 
the  present,  the  MM-GTD  method  is  applicable  under  more  general  situ- 
ations than  those  presently  handled  by  the  GTD,  e.g.,  it  is  applicable 
to  variable  curvature  cylinders  and  variable  surface  impedance  cases; 
furthermore  it  is  applicable  to  the  cases  in  which  the  surface  impedance 
does  not  support  surface  wave,  or  Elliot  type  modes.  Some  refinements  are 
also  required  in  the  W-GTD  methods  as  pointed  out  in  [8];  these  refine- 
ments will  also  be  studied  in  addition  to  the  continuing  GTD  study  of 
the  antenna  problems  of  Figures  1 and  2.  In  conclusion,  the  GTD  method 
does  show  promise  of  being  a useful  method  in  analyzing  the  antenna 
problems  of  Figures  1 and  2. 
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Figure  9.  Radiation  pattern  of  a magnetic  line  source  on  a perfect ly-conau 
planar  surface  which  is  covered  with  an  impedance  surface  patch. 
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Figure  13.  Radiation  pattern  of  a magnetic  line  source  on  a 
perfectly-conducting  circular  cylinder  which  is 
covered  with  an  impedance  surface  patch. 


Mqure  14.  Radiation  pattern  of  a magnetic  line  source  on  a 
perfectly-conducting  circular  cylinder  which  is 
covered  with  an  impedance  surface  patch. 
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Figure  17.  Radiation  pattern  of  a magnetic  line  dipole  source  on  a perfectly-conducting 
planar  surface  wftich  is  covered  with  an  impedance  surface  patch. 


Radiation  pattern  of  a magnetic  line  dipole  source  on  a pertectiy-c 
circular  cylinder  which  is  covered  with  an  impedance  surface  patch. 
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Figure  20.  Radiation  pattern  of  a magnetic  line  dipole  source 
on  a perfectly-conducting  circular  cylinder  which 
is  covered  with  an  impedance  surface  patch. 
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APPENDIX  I 

AN  ANALYSIS  OF  THE  CANONICAL  PROBLEM  OF  THE  RADIATION  FROM 
A LINE  SOURCE  ON  AN  INFINITE  PLANAR  IMPEDANCE  SURFACE 


The  geometrical  configuration  of  the  problem  is  illustrated  in 
Figure  4.  A magnetic  line  source  at  0 excites  this  configuration  of 
interest.  This  line  source  is  designated  by  R,  and  it  represents 
a z-directed  quantized,  equivalent  magnetic  line  source  for  the  axial 
slot  (or  the  TE^)  case;  whereas,  it  represents  an  x-directed  quantized, 
equivalent  magnetic  line  dipole  source  for  the  circumferential  slot 
(or  TM^)  case.  The  TE^  case  is  treated  first;  it  is  followed  by  a 
similar  analysis  for  the  TM^  case. 

TE^  or  the  axial  slot  case 

The  magnetic  line  source  M is  specifically  given  by 


M = z M 6(x)6(y) 


(A-1) 


where  the  strength  M is  assumed  known.  The  magnetic  field  radiated  by 
this  source  in  the  presence  of  the  infinite  planar  impedance  surface 
is  entirely  z-directed  (hence  the  fields  are  TEz  type).  Let  this 
magnetic  field  be  denoted  by  hz.  hz  satisfies  the  reduced,  inhomogeneous 
wave  equation,  and  the  following  boundary  conditions. 


(V^  + k^)  h^  = jkY^  M 6(x)6(y-h] 


h-*-0  and  y^O;  lx|<“. 


(A-;?) 


Tv'  " 7“  ^ 


|xl<“  and  Zg  is  the 

surface  impedance  value.  (A-3) 


One  notes  that  Yq  = free  space  admittance  and  Zq  = Yq"  . is  the  two 
dimensional  Laplacian  operator.  Also,  satisfies  the  radiation 
condition  as  p-«°.  The  solution  to  a similar  problem  is  discussed 
elsewhere  [9,10];  however,  a method  of  solution  based  on  Green's  functions 
is  presented  here  for  the  sake  of  completeness.  Only  the  significant 
steps  in  the  analysis  will  be  indicated  for  the  sake  of  brevity.  Let  g^ 
denote  a Green's  function  which  satisfies  the  following  equations: 


(’t  k‘^)g^(p'lp)  = - 6(x'-x)6(y' -y) ; (y,y')^0,  lx,x'l<“. 


37 


r 


where  is  the  2-D  Laplacian  in  the  primed  coordinate  system  (x',y‘), 
^ = XX  + yy;  whereas  = x'x  + y'y. 


aty-o,  |x' 


(A-5) 


Also,  9m  satisfies  the  radiation  condition.  Utili_ziri£  the  two-dimensional 
Green's  theorem  and  the  reciprocity  condition  9fn(p'lp)  = gm(PlP')* 
obtains 


h,(x,y)  = - jkY  M g(x,ylO,h);  h-^0. 


(A-6) 


Using  the  procedure  for  constructing  higher  dimensional  Green's  function 
as  in  [10],  one  obtains 


- A • cj?  g ( X j X ' ) 
2ttj  X •’’X 


g (y ly ' ) dA 


(A-7a) 


where  the  one  dimensional  Green's  functions  g^x  g^^y  satisfy  the 
following  differential  equations  and  boundary  conditions: 


^7  " \)9mx  ^ VV  = ■ (A-7b;A-7c) 


g and  g satisfy  the 

^mx  ^my  a . ^s 

radiation  condition  at  p-*^;  ^ g„„  = jk  y-  g at  y=0.  (A-7d;A-7e) 

0 ^ 


3y  my 


In  addition 


A + X = k . 
X y 


(A-7f) 


The  contour  Cx  encloses  only  the  singularities  of  g^ix  the  complex  A^^ 


plane,  g^^^^  and  g^^^  are  found  to  be 


""mx 


+ R,  e 


"my 


(A-7q;A-7h) 
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in  which  y > h 0.  The  contour  runs  in  the  counter-clockwise  sense 
around  the  branch  cut  singularity  of  along  ReX^  ^ 0.  Introducing 
the  transformation  = X , one  obtains 

X X 


9„(x.yio,t,)  = jU-  J 


dk 


■J£ 


i ^ + R.  e ^ 

h 

y 


j'^(y+h) 


-jkjxl 


(A-8) 


in  which  x'=0,  y'=h,  and  y^h-K).  Also, 


►X,,  - k Z /Z 
R = -JL 

" /r~  + k Z /Z 

y s'  0 


^ ^s/^0 


(A-9) 


The  contour  of  integration  in  Equation  (A-8)  is  illustrated  in  Figure 
A-1.  'he  branch  cuts  in  Figure  A-1  are  chosen  such  that  Im  /xy<0.  It 
is  convenient  for  subsequent  evaluation  of  Equation  (A-8)  to  introduce 
further  transformations 


Figure  A-1.  Contour  of  integration  in  the  complex  k^  plane. 


X = p cos(fi  ; y = p sini<)  , 


and 


k^  = k cos4 


-‘'x 


+ k sini 


(A-10) 


(A-11) 
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The  Green's  functiun  Oni  tquation  (A-U)  is  now  (?xprHSscrl  in  the 
complex  C plane  as 


/ 1 ^ 1 r ( -jkp,cos(#-4) 

g^(x.ylo.h)  = - ^ j d5  e +R^(Oe 


-jkP2COs(^-4)j  (A-12) 


where 

y + h 5 sin<(i  ; and  the  corresponding  lx|  h cosik.  (A-13a;A-13b) 
2 2 

One  notes  that  p as  h -*•  0,  Also, 

2 


sinC  - 1 II 

= rrn— 


(A-14) 


The  contour  of  integration  c is  indicated  in  Figure  A-2.  The  integral 
in  Equation  (A-12)  may  be  evaluated  for  large  kp-|  via  the  method  of 

2 

steepest  descent.  Thus,  the  contour  c is  deformed  into  the  steepest 
descent  path  SDP  through  the  saddle  point. 


♦c,  ^Cz 


SURFACE  WAVE  POLES 


Kk-+/S 

K«.-)8 


sop;  \ original  contour  c 

(MAPPED  FROM  K.  TO  ( PLANE) 


Figure  A-2.  Contour  of  integration  in  the  C plane. 
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(A-15) 


5 = = ({,  (at  the  saddle  point);  0<5<n/2 

if  0<^<n/2. 

The  case  Q<^<v/2  is  treated  here.  (A  separate  analysis  for  is 

not  required  since  the  results  for  this  case  may  be  obtained  Trom  those 
for  0<(j)<7r/2  by  invoking  symmetry.)  In  deforming  the  contour  c into  SDP, 
the  poles  of  Rh(C)  may  be  crossed,  hence,  their  residue  contribution  must 
be  taken  into  account.  The  poles  of  Rh(^)  occur  at 

sin  = - Z /Z  ; such  poles  certainly  exist  if  Zs  is  chosen 
^ to  be  jX|_  where  Xi_>0  (i.e.,  Zg  is  an 

inductive  reactance).  (A-16) 

= sin'^  {-j  X^^/Z^)  at  the  pole  for  x>0  (A-17a) 

and 

5p^  = TT  + sin'^  (j  X|^/Z^)  at  the  pole  C=5p2  i"or  x<0.  (A-17b) 


Let 


k cos  Cp^  = 8 ; k cos  ^p2  = “ ^ » 


(A-18a;A-18b) 


where 


8 = k^l  + (X^/Zq)  . with  Z^  = JXl(Xl>0). 


(A-19) 


Finally,  a saddle  point  evaluation  of  Equation  (A-12)  yields  the  following 
for  pole  not  close  to  the  saddle  point. 


g - i 


2 sin^  ^ a 


-ay-jelx  j 


I8irk  si  n4i  + Z^/Z  /- 

S'  0 rp 


- J T e 


■]-')>) 


for  0<4><it/2  and  h=0. 


2 2 2 

a H 8 - k , (e^k),  and  6 is  defined  in  Equation  (A-19). 


(A-20) 

(A-21) 


\ 


A 


41 


Also.U  in  Equation  (A-20)  denotes  the  Heaviside  step  function  which  is 
unity  for  whereas  it  is  zero  for  The  angle  <t>c1  is 

illustrated  in  Figure  A-2.  The  first  term  on  the  RHS  of  Equation  (A-20) 
is  the  so-called  space  wave  contribution  which  vanishes  along  the 
impedance  boundary  (at  ^=0,  it);  whereas,  the  second  term  is  the  surface 
wave  mode  field  which  exists  when  = j Xl(Xl>0).  One  may  now  find  h^ 
in  terms  of  Equation  (A-6)  and  (A-20)  to  be 


TM^  or  the  circumferential  slot  case 


The  magnetic  line  source  W at  0 in  this  case  is  a magnetic  line 
dipole  given  by 

M = X M 6(x)6(y)  , (A-28) 


where  the  amplitude  M is  a known  quantity.  The  electric  field  radiated 
by  this  source  is  entirely  z-directed  (hence,  the  fields  are  TM2  type). 

Let  th-’s  electric  field  be  denoted  by  e£.  Then 

+ k^)f  = - M 6(x)6(y-h)  ; y^O;  lxl<»  and  h-K)  (A-29) 

where 

; e^  z = - Vx(xf)  (A-30) 


and 


xf  is  the  electric  vector  potential. 


The  impedance  boundary  condition  bz  z = y x Zs  at  y=0  where  TTtan 
is  the  tangential  magnetic  field  at  y=0  corresponding  to  e^  implies  that 


aty=0,  |x|<-  . 

3y  S 

Let  ge  be  a Green's  function  which  satisfies 

+ k^)gg(F'lD  = - 6(x'-x)6(y'-y)  ; ylO;|xl<» 

y'^0; |x' |<* 


at  y' =0;  jx' i<*  . 


^0  0 

The  above  Equation  (A-33)  could  be  simplified  to  t— r = jk  7-  g„ 

3y  ^5  e 

3 9e 

by  integrating  with  respect  to  y'  and  requiring  that  g and  r 
the  radiation  condition.  One  may  show  that  ^ 


(A-:;i) 


(A-32) 


(A-33) 
at  y'=0 


satisfy 
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f = M gg(plp')  = M gg(x,ylo,h),  h-K). 


(A-34) 


In  obtaining  Equation  (A-34),  use  is  made  of  the  fact  that  ez,  and  f satisfy 
the  radiation  condition;  also,  gg  is  chosen  to  satisfy  the  radiation 
condition.  Following  the  procedure  outlined  for  the  construction  of 
g^  in  the  TE2  case,  one  obtains 


“+je 


g^,(x,yio,h)  = ^ J 


dk^  e 


-jkjxl 


-»-je 


g-jY(y-h)^.[^  g-jY(y+h) 


y>h 

h-Hl, 


(A-35) 


where 


Y 


k -k 

X 


Y-k  Z„/2, 

h'h 


( A-36a;A-36b) 


The  contour  of  integration  in  Equation  (A-35)  is  the  same  as  in  Figure 
A-1 . Transforming  the  above  integral  into  the  contour  integral  over  c 
in  the  complex  C-plane  of  Figure  A-2  via  Equations  (A-10)  and  (A-11),  one 
obtains 


e^(x,y)  = - M ^ 9g(x,y |o,h) 


M ( 
4tt  j 


d5  ksin^ 


-jkp^cos(«ji-5) 


R,e 


-jkppCos((>i-5) 


(A- 37) 


where  p-j  2 defined  in  Equations  (A-1 3a;A-l 3b) , and 


sin5  - Z /Z 

= ■ sWVlV 

0'  s 


(A- 38) 


One  notes  that  Rs(5)  can  certainly  have  a pole  when  Zj  = - j Xc  (Xg^’O); 
i.e.,  when  Zg  is  a capacitive  reactance.  Then,  evaluating  Equation  (A-37) 
by  the  saddle  point  method  while  taken  into  consideration  the  residue 
at  the  pole  of  Rs(5)  when  deforming  c into  SDP  yields  the  following 
result  for  pole  not  close  to  the  saddle  point,  and  for  y-Kl  (i.e.,  p-|  2^9) 
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- jkM  ^ 0i^iV2. 


(A-39) 


The  results  for  -n^^/2  are  directly  obtained  from  the  results  for 
0<(}i<n/2  above  by  Tnvoking  syimetry.  One  notes  that  a and  6 in  Equation 
(S’-TQ)  are  defined  as 

+ k^  (B>k);  a = kZ^/X^  = -jk  Z^/Z^  (with  Z5=-jX^) 

(A-40a;A-40b) 

in  which  Zs  = - jXc(Xg.>0)  implies  capacitive  reactance  as  indicated  earlier. 
One  may  rewrite  Equation  (A-39)  as 


e,(P)  u"^  + u^ 


(A-41) 


where 


= c M R(^) 

*0 


space  wave 


(A-42) 


u^  = c M e j^i^l  = surface  wave  mode  field. 
y=0  ° 


(A-43) 


If  Cq  i;  defined  as 


'0  ' • « 


(A-44) 


R(l)  = = L2  sin^)(jg/k) 

sin(|)  + sin4i  + ja/k 


(A-45) 


I SW  _ /q'  Vl  “ 

L = ''SiTjk  — 


(A-4b) 


in  which  a and  b are  as  defined  in  Equations  (A-40a;A-40b) , with 
Zs  = - 
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APPENDIX  II 

AN  ANALYSIS  OF  THE  CANONICAL  PROBLEM  OF  SURFACE  WAVE 
DIFFRACTION  BY  A PLANAR,  TWO  PART  IMPEDANCE  SURFACE 


The  geometrical  configuration  of  this  canonical  problem  is  il- 
lustrated in  Figure  5.  The  incident  field,  u^  is  a bound  surface  wave 
mode  which  propagates  along  the  surface  impedance  boundary  (y=0,  x<0); 
and,  it  produces  the  reflected  surface  wave,  and  diffracted  fields 
and  u^,  respectively  upon  striking  the  discontinuity  in  surface 
impedance  at  (x=0;  y=0).  One  requires  that  the  surface  impedance  be 
purely  inductive  in  the  TE^  case;  whereas,  it  be  purely  capacitive  in 
the  TMz  case  (to  support  the  surface  wave  mode).  In  the  TE2  case,  the 
magnetic  field  is  entirely  z-directed;  whereas,  in  the  TM2  case,  the 
electric  field  is  entirely  z-directed.  Let 


f 'I 

r \ 

Y ; 

rsw  _ 
u = ^ 

z 

* ; for  < 

TE 

J 

r^SW 

1 ^ J 

™z. 

(A-47a;A-47b;A-47c) 


where  and  refer  to  the  z-directed  electric  and  magnetic  field 
intensities  in  the  TM^  and  TE7  cases,  respectively;  also,  the  superscripts 
i,  rsw  and  d in  (A-47a;A-47b;A-47c)  stand  for  incident,  reflected  and 


diffracted  wave  components,  respectively. 


tutes  a known  excitation,  and  the  fields 


The  incident  field 
-d 


and 


consti- 


can  be  found  via 


the  Wiener-Hopf  technique  for  solving  the  two  part  boundary  value  problem 
of  Figure  5.  The  TE^  solution  will  be  briefly  outlined  first;  and  a 
similar  outline  will  follow  for  the  TM^  case. 


TE^  case 


The  excitation  has  the  form 


h’  = e"^'^  e^^^^  ; |xj^“ 


in  which  it  is  assumed  that  exists  even  for  x>0. 


(A-48) 


46 


6 = k 


1 - 


2 2 


(A-49a;A-49b) 


as  in  Fquations  (A-19)  and  (A-21)  for  the  TE^  case  with  Zj  = jXL(XL>0) 
corresponding  to  inductive  reactance.  In  the  following  analysis  an 
e ^ time  convention  will  be  employed  for  the  sake  of  convenience  in 
using  the  Wiener  Hopf  notation;  howQver,  a complex  conjugate  of  the 
final^esults  will  yield^ack  the  e^"^  time  convention.  The  quantity 
i = r-1  instead  of  j = > -1  will  be  employed  to  distinguish  between  the 
two  time  conventions.  Let  H|  denote  the  scattered  field,  then  the  total 
field,  + H^.  The  total  field  satisfies 


(V^  + k^)  = 0; 


for  y > 0, 


jp  " 7~  -X  ~ X < 0. 

Z 


(A-50) 


(A-bl) 


Equation  (A-51 ) implies  that  ^ = - ik  for  y = 0,  x < 0,  and 


3H 
2 


= 0 


for  y = 0,  and  X > 0. 


(A-52) 


Equation  (A-52)  also  implies  that 


3H' 


3H 


“3y  " - -57  l^or  y = 0,  X > 0. 

satisfies  the  radiation  condition  for  e time  dependence. 
One  may  next  define  the  following  Fourier  transforms. 


(A-53) 


h^  = 1 


*/ 


1SX 


e'""  dx  = h;  + h^ 


uS  „isx  . rs  1 
H e dx  ; h = — i— 


(A-54a) 


u 

/ 


e^^^  dx 


Ob 

i _ 1 r Ljl  i 

* ■ 'L  ' 

K 2tt^ 


sx  , ie 
dx  = — 


-ay 


* .’^(s+d) 


, via  (A-48). 


(A-54b;A-54c) 


(A-55) 
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h|  and  h1  individually  satisfy  the  wave  equation  in  Equation  (A-50). 
Fourier  transforminq  Equation  (A-50)  for  instead  of  H^,  one  obtaines 


(A-56) 


in  which  Im  /x>0'.  The  solution  to  Equation  (A-56)  in  accordance  with 
Equation  (A-53)  is 

h^  = c(s)  , y 1 0.  (A-57) 

^ s s 

The  inverse  Fourier  transform  of  h leads  to  Thus, 


00 


(A-58) 


Incorporating  Equations  (A-48)  and  (A-58)  in  H2  = which  appears 

in  the  boundary  conditions  (A-51)  and  (A-52)  yields 


00 

j L"^ (s)^(s)e'^^^ds  =0  , X < 0 , (A-59) 

^ 00 


l>{s)e 


-isx 


ds  = - ii'ZTT  a e 


iex 


X > 0 


(A-60) 


where 


!f(s)  s y'x  c(s) 


and 


L(s)  = 


✓ A 

/T  - i 


la 


Nk^-s' 


I 2 2 

k -s  - ia 


1 


(A-61) 


(A-62) 


The  above  set  of  dual  integral  equations  for  (|)(s)  Tn  Equation  (A-59) 
and  Equation  (A-60)  may  be  solved  via  the  method  of  factorization  [12]. 
Without  going  through  the  details,  one  obtains 
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defines  the  Wiener  Hopf  factors  L (s)  and  L_(s)  for  the  L(s)  of  Equation 
(A-62).  L+(s)  is  analytic  for  Im  s>0  and  L_(s)  is  analytic  for  Im  s<0. 
Incorporating  Equation  (A-63)  in  Equation  (A-57)  and  employing  Equation 
(A-58)  yields 


^ 

Z Ztt 


L_(s) 

L T-e) 


1 


(s+e)  ,k^-s^ 


y e-'isx 


ds 


(A-65) 


The  residue  of  the  pole  at  s=6  in  the  integrand  of  Equation  (A-65) 
furnishes  the  reflected  surface  wave  for  x<0,  whereas  the  pole  at  s=-e 
furnishes  the  negative  of  the  incident  field  for  x>0  (see  Equation  (A-43) 
where  it  was  assumed  that  exists  for  lx|<“),  and  therefore  cancels 
the  incident  field  for  x>0.  The  reflected  surface  wave  field  is 


H 


where 


rsw 


R* 


X < 0 


(A-66) 


2e[L^(0)]' 


Ru  = lim  (s-b)  L(s)  = ^ 
" s-^B  ^ 


(A-67) 


Clearly,  is  the  reflection  coefficient  at  (x=0,  y=0)  associated  with 
the  reflected  surface  wave. 

The  diffracted  field,  H7  may  be  found  from  an  asymptotic  evaluation 
of  the  integral  in  Equation  (A-65)  via  the  method  of  steepest  descent. 

Jo  this  end,  one  introduces  the  transformations  s=kcosw,  ds=-ksinwdw, 

* 2 2 ^ * 
ijk  -s  =ksinw,  x=pcoS(|i  and  y=psin(>.  Thus, 


= 


h f 


L (kcosw)  / . X 

- ^-ikpcos(w+4i)  . 

L_(kcosWp)(kcosw+B) 
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(A-68) 


The  contour  in  Equation  (A-68)  is  indicated  in  Figure  A-3.  The 
contour  c^^  may  be  defonned  into  the  steepest  descent  path  (one  notes 
that  the  pole  at  s=+6  which  may  be  crossed  in  the  path  deformation  yields 
of  Equation  (A-66)  and  hence  this  pole  crossing  will  be  ignored 
since  it's  contribution  has  already  been  evaluated  in  Equation  (A-66)), 
Without  any  details  for  brevity,  the  resulting  saddle  point  approximation 
along  the  steepest  descent  path  yields 


, a j"2Tt 
^ kp 


(kP-  J) 


L(-kcos^) 

L^.(-kcos<ti)L^(kc( 


;e-kcosi(>J 


The  quantity  Wp  in  Equations  (A-68)  and  (A-69)  is  related  to  8 via 

8 = kcosw  . 

P 

One  may  rewrite  Equation  (A-69)  as 


'x-  D*(^)  ^ 

/p 

where 


D*(«)  = - 


L(-kcos<|)) 


L^(-kcos(()] 


; B-kcoS(() ' 


The  explicit  form  of  the  factors  L±(s)  may  be  obtained  conveniently 
from  the  procedures  outlined  in  [13]  or  [14],  In  the  present  instance, 
it  has  been  found  that  the  procedure  given  by  Weinstein  [14]  leads  to 
simple  and  useful  expressions  for  Lt(s).  Without  going  through  the 
details  of  obtaining  L±(s)  via  the  method  in  [14],  one  may  directly 
write  the  expressions  for  Rh  and  D(ifi)  corresponding  to  the  surface  wave 
reflection  and  diffraction  coefficients  in  the  eJ^^  time  convention  in 
terms  of  the  results  for  L±(s)  as 


si  nh  u 


1-8/k 

~8A” 


(A-73) 
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.7T 

Jj  

D(^)  = S . 2coshs(Hcos»)  . 

cosh^-coS(^  ( 1 +coshC ) { coshs+cosifi ) 


with 


sinh 


0^ 

k 


♦-jC 

udu  ^ 1 r udu 

sTnfiir  7n  J sTruT 


(A-74) 


(A-75) 


It  is  noted  that  Rh  and  D in  Equations  (A-73)  and  (A-74)  for  the  e^*^^ 
time  dependence  are  the  complex  conjugates  of  and  D*  in  Equation  (A-67) 
and  Equation  (A-71)  for  the  time  dependence,  respectively. 


TM^  case 


The  excitation  E^  in  this  case  is  given  by 

ri  _ „-ay+i3x  I I 

E^  = e ^ ; |x|<» 

in  which  it  is  assumed  that  E^  exists  even  for  x > 0. 


2 ? 2 

a = - k (8>  k)  , 


(A- 76) 


(A-77a;A-77b) 


as  in  Equations  (A-40a;A-40b)  for  the  TM^  ca'.e.with  Zs=-jX(-(Xc>0)  corre- 
sponding to  capacitive  reactance.  Again  an  e"^“^  time  dependence  will 
be  employed  in  the  analysis  to  follow;  however,  the  final  results  will 
be  presented  for  the  eJ*^^  time  dependence.  Let  E^  and  e!. denote  the 
total  and  the  scattered  fields,  respectively,  then  E7  = e1  + Ef.  The 
total  field  E^  satisfies 


(V^  + k^)  E^  = 0 


; for  y ^ 0,  |x| 

; for  y = 0 and  x 
; for  y = 0 and  x 


radiation  condition 


< 00 

9 

< 0 , 

> 0. 
for 


(A- 78) 

(A-79) 

(A-80) 

time  dependence. (A-81 ) 


E^  satisfies  the 
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One  defines  the  Fourier  transforms 


^s 

e = 


— r e^^^  dx  = e^  + as  in  (A-54a) 
2n  J ^ 


yZn 


and 


QP 

— r E^ 


isx  _ 1 e 
e dx  = 


-ay 


/2tT  ^ /2V(s+e) 


, via  (A-76) 


e|  and  E2  individually  satisfy  the  wave  equation  in  (A-78).  Thus, 
Fourier  transforming  Equation  (A-78)  for  e|  instead  of  E^  yields 


= ' k^-s^  = i 


i)|s^-k^ 


in  which  Im/x  > 0.  The  solution  to  Equation  (A-84)  subject  to 
Equation  (A-81 ) is 


i/xy 


e = A(s)  e ^ 1 y 1 0. 


The  inverse  Fourier  transform  of  Equation  (A-85)  leads  to  E' 


e!  = ^ 


/2V 


00 

-J 


e^  e-^^"  ds 


One  may  incorporate  Equations  (A-76)  and  (A-86)  in  the  boundary  con^ 
ditions  of  Equations  (A-79)  and  (A-80)  to  obtair; 


00 

/ 


G'''(s)A(s)  e"''^^  ds  = 0,  x < 0 


/ 


A(s)  ds  = - e’“^  , x > 0. 


where 


G"^(s)  = /x  -ia  = fX  L"^(s)  , (see  (A-62)  for  L(s)). 


b2 


(A-82) 


(A-83) 


(A-84) 


(A-85) 


(A-86) 


(A-E7) 


(A-88) 


(A-89) 


Equations  (A-87)  and  (A-88)  are  a set  of  dual  integral  equations  for 
A(s)  which  may  be  solved  via  the  method  of  factorization  as  in  TEz  case 
treated  earlier.  Without  presenting  any  details  for  the  sake  of  brevity, 
the  result  for  A(s)  is  given  by 


1 2__(s) 

Ts^  / (-eT 


(A-90) 


where.  G(s)  = G+(s)G.(s)  in  which,  G+(s)  and  G_(s)  represent  the  Wiener- 
Hofp  factors  of  G(s).  The  A in  Equation  (A-89)  may  be  factorized  by 
inspection,  and  the  factorization  of  L(s)  is  known  from  the  TEz  case; 
thus  the  factorization  of  G(s)  is  directly  obtained.  One  may  now  write 
Equation  (A-86)  as: 


-1 

2n 


G_(s) 

G~(-e) 


TsW 


y j-'s* 


ds 


(A-91) 


As  in  the  TEz  case,  the  pole  of  the  integrand  in  Equation  (A-91)  at 
s=-g  yields  a residue  contribution  only  for  x>0  which  exactly  cancels 
the  incident  field  for  x>0.  On  the  other  hand,  the  pole  at  s=6  yields 
the  reflected  surface  wave  E^^'^  for  x<0. 


^ ^-ay  g-iex 


X < 0 


where 


R|  = 


26[G^(B)] 


7 


R = lim  (s-e)  G(s)  = 

® s-^e 


ia 

6 • 


(A-92) 


(A-93) 


R£  represents  the  reflection  coefficient  at  (x=0,  y=0)  for  the  reflected 
surface  wave  field. 

The  diffracted  field  Ez  may  be  obtained  by  asymptotically  approxi- 
mating Equation  (A-91)  via  the  method  of  steepest  descent  as  done 
earlier  for  in  the  TE2  problem.  Thus, 


^ D*(^) 


+ikp 


where 


D* 


G(-kcos<|i)  sini^  . 

G^T-kcoS(Ji  )G^(  KCosWp)  ( e-kcos4i ) ’ 


e=kcosWp. 


(A-94) 


(A-95) 
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As  in  the  TE^  case,  the  surface  wave  reflection  and  diffraction 
coefficients  Rf  and  D may  be  explicitly  written  down  for  the  eJ“^  time 
convention  as 


,(3/k)^  - 1 


TT  J s^i  nViu 


(A-96 


D(a)  = i ^ Sint  ] J ^coshS  0 

V‘^°sh?-cost/  'IcoshC+cost 


j_  / udu  + I_  [ udu 

2tt  / sTnTur  2tt  / sinu 


(A-97 


where 


5 = sinh 


._u-l  o 


(A-98 


R£  and  D for  the  e'^^^  time  convention  are  the  complex  conjugates  of 
R|^  and  D*  in  Equations  (A-93)  and  (A-95)  for  the  case,  respectively 


Im 


Re  w 


Figure  A-3.  Contour  of  integration  in  the  complex  w-plane. 


APPENDIX  III 


SURFACE  FIELDS  IN  THE  CANONICAL  PROBLEM  OF  A LINE 
SOURCE  ON  A CIRCULAR  CYLINDER  COMPLETELY  COVERED 
WITH  AN  IMPEDANCE  SURFACE 


The  geometrical  configuration  of  the  problem  is  illustrated  in 
Figure  7.  The  line  source  is  a magnetic  current  fT  which  is  either  a 
z-directed  magnetic  line  source,  or  a ^-directed  magnetic  line  dipole; 
as  before  the  former  represents  the  source  for  the  TE^  case;  whereas, 
the  latter  represents  the  source  for  the  TM2  case.  The  TE,  case  is 
treated  first,  and  it  is  followed  by  the  treatment  of  the  TM2  case. 

TE^  or  the  axial  slot  case 


The  magnetic  current  M at  Q'  in  this  case  is 
M = z M 


6(p-a)i5  (») 

p 


(A-99) 


and  it  generates  only  a z-component  of  the  magnetic  field  which  is 
denoted  by  h2. 

(V?  + k^)h  = ikY  M ; p > a 

''  ^ ° P 0 ~ ^ < 2tt  . 


and 


3h 

3p 


- - jk  ^ 

0 


0 at  p = a. 


(A-lOO) 


(A-101) 


Also,  hz  satisfies  the  radiation  condition.  Using  the  method  of 
Green's  functions,  one  may  obtain  the  following  expression  for  h2; 


h^lP,^)  = - jkY^  M G^(pIp') 


(A-102) 


p'=a 


where  Gn,  (pIp')  is  a Green's  function  which  satisfies  the  radiation 
condition,  and 


(v;^  + k^)  G (^'  1^")  = - 
' t ' m'  ' ' 0 


in  which  p is  the  position  vector  of  (p,4>);  and  likewise  p'  is  the  position 
vector  of  (p  ' ,<j> ' ) . Also, 
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= 0 at  p = a 


(A-104) 


; Z = value  of  the  surface  inpedance  which 
^ completely  covers  the  cylinder; 

Zq  = free  space  impedance. 


(A-105) 


The  solution  to  may  be  constructed  via  the  methods  presented  in  [10]; 
thus 


37" 


C0Sv(7t-  I 

si  nvii 


(A-106) 


where  H'‘^'(kp)  is  the  cylindrical  Hankel  function  of  the  second  kind, 
and  of  order  v.  One  may  write 


cosy(TT- 1 
si  nvir 


” j g-jv(2ffA) 

1=0  U -* 


(A-107) 


Physically,  the  above  series  corresponds  to  multiple  encirclements  of  tne 
field  around  the  cylinder  in  the  azimuthal  propagation  representation 
for  Green's  function  Gfp  presented  in  Equation  (A-106).  Since  one  is 
interested  in  applying  the  results  of  the  present  analysis  to  perfectly 
conducting  cylinders  with  a surface  impedance  cover  of  finite  extent,  the 
effects  of  these  multiple  encirclements  will  be  neglected.  Thus,  for 
the  present  purposes  one  writes 

Gj^(pjr')  = Sp^(p'|^')  + multiply  encircling  wave  contributions  (A-108) 


from  which  it  follows  that 

1 dv  ^ 

«'=0  ([4--jk4]Hf^(kp)[ 


-jkA]H^)(kp) 


(A-109) 


and  the  field  h^  which  corresponds  to  this  non-encircling  part  is  h2< 
Thus, 


56 


h,(p.«)  = - jkY„  M S (p.^la.O) 


(A-no) 


represents  the  non-encircling  part  of  hz  in  terms  of  in  which 
Gm(p,*;a,0)  is  given  explicitly  in  Equation  (A-109).  Let 


h^(p,*)  = + h^ 


(A-111) 


where  the  h2  terms  corresponds  to  the  clockwise  circumferential  propagation 


cue  no  te 

.-jvlOl  i. 


Equation  (A-109),  whereas  h^  corresponds. to  the  counter 
ferentially  propagating  term  e“Jv(2ir- l(j) ] ) -jn  Equation 


term  e 

clockwise  circumferentially  propagating 
(A-109),  Defining 


2it-  |iji  I 


m\|i 


± /ka\ 

; m = 


1/3 


(A-112;A-113;A-114) 


and  assuming  that  ka  is  sufficiently  large,  one  may  approximate  and 
hence  n|  on  the  surface  p=a,  within  the  shadow  region  of  the  source  at 
Q'  (p=a‘,‘  (j>=0)  as 


p=a 


,-jka* 


dt 


W2(x)e 


■jcS 


W2(T)+jm  j- 


(A-115) 


in  which 


*^0  = 


'ii?  ^0 


(A-116) 


and 


H(2)(ka)  (t);  ^n(2)(kp)  ^ - ii^w'(x), 

m/iT  mr  TT 


(A-117a;A-117b) 


The  variables  v and  x are  related  through  the  transformation: 


_ V - ka 

X = 


m 


(A-117c) 


W2(x),  and  wMx)  are  the  Fock  type  Airy  functions  defined  in  [5] 

(also  see  [6]).  If  the  surface  impedance  Zg  is  inductively  reactive, 
then  the  large  curved  impedance  surface  will  support  an  Elliot  type 
mode  mentioned  in  Sections  I and  II.  The  field  of  this  Elliot  mode  is 
obtained  by  evaluating  the  residue  of  the  integrand  in  Equation  (A-115) 
at  x=xq  which  locates  the  Elliot  mode  pole  of 


L_ 
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1 

z 

W2(t)+jin  ^ W2(t) 


Thus, 


(A-118) 


While  several  papers  have  been  written  on  the  solution  for  tq  in 
Equation  (A-118),  the  one  which  appears  to  have  dealt  with  this  problem 
in  great  detail  is  due  to  Logan  [5];  and,  Logan's  results  are  used  in  the 
present  work.  The  value  of  tq  is  approximately  given  by  [5] 


X,  ^ (q‘ 


+ + 


J 

8q 


32q' 


_n 

32q 


) . J 2,2  e-f/S?' 


q»2,  (A-119) 


where 


In  addition  to  the  root  t=Tq  corresponding  to  the  Elliot  mode,  there 
are  other  roots  of  the  equation  in  (A-118)  which  correspond  to  the  set 
of  Watson  modes  which  are  characterized  by  a significant  exponential 
decay  along  their  propagation  paths.  One  is  primarily  interested  in  the 
Elliot  mode  field  in  this  study  as  mentioned  in  Section  II-B,  hence  the 
Watson  type  modes  will  not  be  discussed.  Evaluatino  the  residue  at 
T=TQ  yields  the  Elliot  mode  field  corresponding  to  h*  as: 


h^la,^.) 


kV  M 

0 e 


which  may  be  rewritten  as 
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(A-121) 


(A-122) 
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(A-123) 
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PW 

The  quantities  L and  y are  given  by: 


Z ■ 
s 

T 

0 


and 


(A-124) 


in  which  m and  tq  were  previously  defined  in  Equation  (A-114)  and 
Equation  (A-119),  respectively,  Cq  is  defined  in  Equation  (A-116). 

The  space  wave,  or  the  field  directly  radiated  by  the  source  may 
be  found  by  a stationary  phase  evaluation  of  the  integral  in  Equation 
(A-106)  for  the  lit  region, after  replacing  cosv[n- l(j>  1 J in  Equation 
(A-106)  by  e"J^^cosv(J>+jeJ^  1'^’ I sinvTi,  and  by  including  only  the  term 
eJ'^l'Hsinvn  ip  this  asymptotic,  stationary  phase  evaluation.  The 
other  term  e"J'^'^cosv(<i  gives  rise  to  diffraction  contributions  arising 
from  the  multiply  encircling  circumferentially  propagating  modes  which 
are  not  of  interest  in  this  analysis  as  discussed  earlier  in  the 
evaluation  of  Gm  or  h^  in  Equation  (A-109)  and  Equation  (A-110).  This 
stationary  phase  evaluation  leads  to  the  geometrical  optics  approximation 
for  the  field  directly  radiated  by  the  slot  as 


c M R(iji) 
lit  region  ° 


/s^ 


s'=Q'P 

Q'  is  at  (a,0) 

P is  at  (A-125) 


where  Cq  is  as  in  Equation  (A-116),  and  R(()()  is  the  same  as  in  Equation 
(A-26)  for  the  planar  case  with  the  exception  that  <(i  in  Equation  (A-125) 
above  is  the  complement  of  the  angle  (ji  in  Equation  (A-26)  due  to  a 
rotation  of  the  coordinate  systems  in  Figures  4 and  6 (or  7). 


TM^  or  the  circumferential  slot  case 


The  magnetic  current  M at  Q'  in  this  case  is 

H = ♦ M . 

p 

which  is  a magnetic  line  dipole  at  Q' , and  it  generates  only  a 
2-component  of  the  electric  field  which  is  denoted  ^ The  field  e^ 
may  be  expressed  in  terms  of  a Green's  function  Ge(p'lp)  which  satisfies 
the  radiation  condition,  the  partial  differential  equation  in  (A-103), 
and  the  boundary  condition  in  Equation  (A-106)  with 
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- e 


(A-132) 


in  which 


(A-133) 


If  the  surface  impedance  Z5  is  purely  capacitive,  then  the  large  curved 
impedance  surface  will  support  an  Elliot  type  mode  mentioned  in 
Sections  I and  II.  The  Elliot  mode  corresponds  to  the  value  of  t=tq 
for  which 


Z 

w^Cto)  + j rn  ^ 


(A-134) 


This  value,  t=tq  corresponds  to  the  location  of  the  pole  of  the  Elliot 
mode  in  the  integrand  of  Equation  (A-132).  The  root  t=Tq  is  obtained 
via  Equation  (A-119)  given  by  Logan  [5]  with 


q = - j 


(A-135) 


Other  roots  of  Equation  (A-134)  corresponding  to  a set  of  Watson  modes 
will  not  be  discussed  due  to  reasons  mentioned  in  Section  II.  Evaluating 
the  residue  at  t=tq  in  Equation  (A-132)  leads  to  the  following  result 
for  the  Elliot  mode  field. 


^ kM  0 e ^ ° 
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(A-136) 


which  may  be  rewritten  as 


-jvt. 


e^(a.^)  = Cq  M e 


; = aij; 


(A-137) 


The  quantities  L and  y are  given  by 


61 


F 


1 

( 


I fiw  /q-  . 0 1 

L = ^ j ^ s • 


'o1 
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Y = k + — T„ 

a 0 


(A-138) 


(A-139) 


The  Cq  in  Equation  (A-137)  is  defined  previously  in  Equation  (A-133), 
and  Tq  is  as  in  Equation  (A-119)  with  q as  in  Equation  (A-135). 

The  evaluation  of  the  field  directly  radiated  by  the  source  in  the 
lit  region  is  not  described  since  it  reduces  to  the  geometrical  optics 
^approximation  as  in  the  TEz  case.  Denoting  this  result  for  ez(p,$)  by 
e^(p,ij>)|  , one  may  write 

I lit  region 


e,(p,4') 


^ c M R((|>) 

lit  region 


1 


Q'P  = s' 

Q'  is  at  (a,0) 

P is  at  (p,4i).  (A-140) 


in  which  Cn  is  as  in  Equation  (A-137),  and  R(^)  is  the  same  as  in 
Equation  (A-45)  for  the  planar  TM2  case  except  that  the  in  Equation 
(A-140)  is  the  complement  of  the  angle  in  Equation  (A-45)  due  to  a 
rotation  of  the  coordinate  systems  in  Figures  4 and  6 (or  7). 
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APPENDIX  IV 

THE  RADIATION  BY  AN  EQUIVALENT  LINE  SOURCE  AT  A DISCONTINUITY 
IN  SURFACE  IMPEDANCE  ON  A PERFECTLY  CONDUCTING  CIRCULAR 
CYLINDER  COVERED  WITH  A SURFACE  IMPEDANCE  PATCH 


\ 

I 


I 


The  geometry  of  this  problem  is  illustrated  in  Figure  8. 
Specifically,  the  problem  to  obtain  an  approximate  but  accurate 
expression  for  the  fields  u^  and  u^  which  are  produced  via  dif- 
fractiojij  of  the^incident  field  u"*  from  Qi.  A similar  s^t  of  diffracted 
fields  U2  and  uX  are  also  produced  via  diffraction  of  u’  from  Q2,  but 
for  the  sake  of  definiteness,  one  may  analyze  only  the  diffraction  from 
Ql  in  this  appendix;  the  results  for  the  diffraction  from  Q2  would  be 
similar  to  those  for  the  diffra;^tion  from  Q] . For  an  electrically  large 
cylinder,  the  diffracted  field  u*]f  on  the  lit  side  of  SBQ]  in  the  deep 
lit  region  is  taken  to  b^  characterized  by  the  same  diffraction  coef- 
ficient as  in  the  field  u°  which  occurs  in  the  problem  of  surface  wave 
diffraction  by  a two  part  impedance  surface  of  Figure  5.  The  diffraction 
coefficient  D is  defined  in  Equation  (A-74)  for  the  TE^  c^se,  and  in 
Equation  (A-97)  for  the  TM7  case,  in  Appendix  II.  Thus,  u?  may  be 
written  as  in  Equation  (20). 


^ u\Qi)  D(6^) 


-jks 

e 


1 


s^  = Q^P  and  P is  on  the  lit  side 
of  SBQ-|  in  the  deep  lit  region. 


(A-141) 


The  field  i^lC^gniay  be  assumed  to  be  produced  by  an  equivalent  magnetic 
line  current  at  Q-j  on  a perfectly  conducting  or  "unperturbed" 
cylinder;  i.e.,  on  the  same  cylinder  as  in  Figure  8,  but  without  the 
surface  impedance  patch.  Onc^  is  found,  it  may  be  systematically 
employed  to  obtain  the  field  Uj  on  the  shadow  side  of  SBQ].  The 
TEz  or  the  axial  slot  case  will  be  analyzed  first;  the  TMz  or  the 
circumferential  slot  case  will  be  analyzed  subsequently. 


TE^  case 


magnetic  field  is  entirely  z-directed  in  this  case,  thus  the 
(P)  in  Equation  (A-141)  may  be  considered  to  be  produced  by  a 
z-directed  magnetic  line  source  at  Q'  on  the  "unperturbed"  cylinder 
as  mentioned  above.  In  this  case,  M®*^  is  given  by 


field  u? 


= z ,,eq 


6(o-a)6(«-$^) 

- 


at 


(A-142) 
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The  z-directed  magnetic  field  produced  by  of  Equation  (A-142)  at  Q' 
on  the  unperturbed  cylinder  may  be  easily  found  to  be  given  by  [6] 

2 c^  ^ ; s-j  = Q^P,  and  P is  on  the  lit  side 

"^^l  of  SBQ^  in  the  deep  lit  region  (A-143) 

where 


and  = 1 


(A-144a;A-144b) 


Comparing  the  RHS  of  Equation  (A-141)  with  Equation  (A-143),  one  obtains 


u'(Q^) 


D(6^)  = 


(A-145) 


in  which  u^(Qt)  and  D(6])  are  given  in  Equation  (15)  (also  see  Equations 
(A-123)  and  (A-124))  and  (A-74),  respectively.  Consequently,  depends 
on  the  aspect  angle  6],  however  this  is  acceptable  since  is  only 
an  equivalent  line  source.  The  value  of  Mf^(6'])  at  6i=0  would  provide 
the  strength  of  the  equivalent  source  which  launches  the  surface  ray 
modes  corresponding  to  u^'^  shown  in  Figure  8.  The  diffraction  from  the 
surface  ray  modes  in  turn  produces  the  field  as  in  Figure  8.  Thus, 

2(l  may  be  directly  obtained  via  the  results  in  [6]  as  indicated  in 
Equation  (22). 
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A. 


P in  the  deep  shadow  region  on  the 

shadow  side  of  SBQi ; A]=QiTi; 

see  Figure  8,  (A-146) 


where 


Lp(Q^)  = - j(jk 


(ka)  Dp, 


(A-147) 


and  Lp,  Dp,  and  ap  respectively  correspond  to  Lp,  Dp,  and  in  [6]. 

In  the  transi tion^region  adjacent  to  the  shadow  boundary  SBQi , one  must 
employ  the  Fock  approximation  in  terms  of  g(*)  as  indicated  in  Equations 
(24)  and  (25). 
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F 


TM^  case 


The  electric  field  is  entirely  z-directed  in  this  case,  thus  the 
field  uj(P)  in  Equation  (A-141)  be  considered  to  be  produced  by  a 
♦o  directed  magnetic  line  dipole  1^*1  at  Q'  on  the  "unperturbed" 
cylinder  as  mentioned  above.  Thus,  in  this  case,  is  given  by; 


= 


‘o 


6(p-a)6((j)-())^) 


at  Q^. 


(A-148) 


The  z-directed  electric  field  produced  by  in  Equation  (A-148)  at 
Q'  on  the  unperturbed  cylinder  is  known  to  be  [6] 


2 c^  ■ ; s^  = Q^P  and  P is  on  the  lit  side 

^1  of  SBQ^  in  the  deep  lit  region,  (A-149) 


where 


jk 


^0  - ■ R ^(Si)  = sin6^ 


(A-150);(A-151) 


One  may  now  find  by  comparing  Equation  (A-149)  with  the  RHS  of 
Equation  (A-141).  Thus, 


,eq  “'(Ql) 

^1  ■ JT^  • 


(A-152) 


in  which  u^(Qi)  and  D(6i)  are  given  in  Equation  (15)  (also  see  Equations 
(A-138)  and  (A-139))  ana  (A-97),  respectively.  Thus,  in  Equation 
(A-152)  is  also  dependent  on  6i  just  as  for  the  TE^  case  discussed 
earlier.  Following  a similar  line  of  reasoning  as  indicated  for  the 
TEz  case,  the  value  of  M^^^at  6]=0  provides  the  strength  of  the  equivalent 
source  which  produces  the  field  in  the  deep  shadow  region  on  the 
shadow  side  of  SBQi  as  shown  in  Figure  8.  One  notes  from  Equation  (A-152) 
that 


= lim 


u'(Ql)  0(6^)  U^Q^) 


30(6^) 


J«i=o 


(A-153) 


It  is  noted  from  Equation  (A-97)  however  that  D(6i)  is  explicitly  pro- 
portional to  sinoj  and  hence  the  value  of  Mp(6-|=0)  is  trivially 
calculated.  Finally,  from  [6], 
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where 
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•^s 


P in  the  deep  shadow  region  on  the 

shadow  side  of  SBQ] ; 4i=QiT]; 

see  Figure  8,  (A-154) 


1/2  / o \ < 

- (jk  f)  (ka)  Dp. 


(A-155) 


Lp,  Dp,  and  cip  respectively  correspond  to  Lp,  Dp  and  ap  in  [6].  In 
transition  region  adjacent  to  SBQ]  one  must  employ  the  Fock  approxi- 


and  Lr 

the  transition  regifjn  — — 

mation  in  terms  of  g(*)  as  indicated  in  Equations  (24)  and  (25). 
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metrh:  system 


BASE  UNITS: 


C^uaniit} 

1 'nil 

SI  Symbol 

Formula 

iength 

metre 

m 

mass 

kilogram 

timp 

sec  ond 

s 

elei  trie,  f urrent 

ampere 

A 

thermiidynamir  tempftratur»* 

kelvm 

K 

amount  of  substance 

mole 

mol 

luminous  intensity 

1 andela 

cd 

SUPPLEMENTARY  UNITS 

pidne  angle 

radian 

rad 

eulid  angle 

sleradian 

sr 

DERIVED  UNITS 

A(  celeration 

metre  per  se<  ond  squartxi 

m/s 

ji  tivily  (of  a radioac  five  soun  ei 

(fisintegration  per  .second 

(difintegration)/s 

angular  acceleration 

radian  per  sec  ond  squared 

rad/s 

angular  velocity 

radian  per  second 

rad/s 

rtftd 

.s(|uare  metre 

m 

di-usitv 

kilogram  per  c uhii  metre 

kg/m 

eledin  lapacitance 

farad 

K 

AsW 

i‘le<  txi(  al  conductance 

sic'mens 

S 

AA/ 

elfi  tru  field  strength 

volt  per  metre 

V'm 

electric  inductaru  e 

henry 

M 

V-s/A 

etc»ctric  potential  difference 

volt 

V 

W/A 

electric  resistance 

ohm 

V/A 

elertromotive  force 

volt 

V 

W/A 

energy 

loute 

1 

N-m 

entropy 

ioulc'  per  ktdvif. 

|/K 

force 

n»?wton 

N 

kg-m/s 

frcH^ueni  y 

hertz 

Hz 

(cycle  )/f 

illuminance 

lux 

lx 

Im/m 

iwnunartce 

candela  per  square*  metre 

cd/m 

luminous  flux 

lumen 

Im 

cd*ar 

magnetic  field  strength 

ampere  per  metre 

A/m 

magnetic  flux 

weber 

Wb 

V-s 

magnetic  flux  densKv 

tesla 

T 

Wfa/m 

magnetomotive  Ion  e 

ampere 

A 

p u\er 

watt 

W 

I's 

f ressuru 

pa.si  al 

Pa 

N/m 

(uantity  of  ele<  trie  ity 

coulomb 

C 

A'S 

(Quantity  of  heat 

imile 

I 

N-m 

radiant  intensity 

e.all  per  sleradian 

W/sr 

specify  heat 

i<iiile  per  kilogram  kelvin 

I'kg-K 

stress 

pasra) 

Pa 

N/m 

thermal  ( onductiviiv 

Walt  per  metre  kels  in 

W/mK 

velo<  ily 

nieire  per  sei  ond 

mil 

viscosity  dynamic 

pascal-second 

Pa-s 

MMosity  kinematif 

square  metre  per  setond 

m/s 

Voltage 

volt 

V 

W/A 

’.filumr 

( ulm  metre 

m 

’AavenumU*r 

rei  iproral  metre 

(wavr)/m 

••.ork 

loiile 

N-m 

SI  PREFIXES 


•iijmi  ( ill  lure  I’rxflx  SI  Syrobol 


1 000  000  oooooo 

) 0 ' ^ 

litra 

T 

1 000  000  000 

gig* 

(; 

1 000  (MH) 

1 1 |S> 

mega 

M 

1 OtMl 

III* 

kilo 

k 

itto 

1 0 ^ 

hecio* 

h 

10 

deka* 

d« 

n 1 

1 0 * 

deci* 

d 

0 ni 

1 0 ' 

i.nnll* 

c 

0 00) 

mini 

m 

0 000  (MM 

mli:iT) 

M 

0 000  (too  001 

nano 

n 

0 (MM)  000  000  (MM 

1(1  ' ^ 

|)l(.0 

|> 

0 (MtO  (too  (mo  000  (HM 

10  ' ^ 

Inmlti 

1 

O (i(  (1  (HMt  000  (MMl  000  (MM 

10  ' * 

alln 

I ■ inii'li  utiiTc 


MISSION 

of 

Rome  Air  Devebpment  Center 


RAVC  plans  and  conducts  rasearch,  exploratory  euid  advanced 
dewlopmant  programs  in  ooaaaand,  control,  and  communications 
fC^J  activities,  and  in  the  areas  of  information  sciences 
and  intelligence.  The  principal  technical  mission  areas 
are  coessunications , electromagnetic  guidance  and  control, 
surveillance  of  ground  and  aerospace  objects,  intelligence 
data  collection  and  handling,  information  system  technology, 
ionospheric  propagation,  solid  state  sciences,  micromeve 
phjrsics  and  electronic  reliability,  smintainability  and 
compatibility. 
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